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Abstract. Over the past few years quadratic Backward Stochastic Differential Equations (BS- 
DEs) have been a popular field of research. However there are only very few examples where 
explicit solutions for these equations are known. In this paper we consider a class of quadratic 
BSDEs involving affine processes and show that their solution can be reduced to solving a system 
of generalized Riccati ordinary differential equations. In other words we introduce a rich and 
flexible class of quadratic BSDEs which are analytically tractable, i.e. explicit up to the solution 
of an ODE. Our results also provide analytically tractable solutions to the problem of utility 
maximization and indifference pricing in multivariate affine stochastic volatility models. This 
generalizes univariate results of Kallsen and Muhle-Karbe [40] and some results in the multivari- 
ate setting of Lcippold and Trojani [42] by establishing the full picture in the multivariate affine 
jump-diffusion setting. In particular we calculate the interesting quantity of the power utility 
indifference value of change of numeraire. Explicit examples in the Heston, Barndorff-Nielsen- 
Shcphard and multivariate Heston setting are calculated. 



1. Introduction 

Since Bismut [9] introduced linear BSDEs in the context of Pontryagin's maximum principle, 
they have been intensively studied. Their popularity stems from the fact that they can be applied 
to many different areas, e.g. in the study of properties of partial differential equations (PDEs), see 
Briand and Confortola [11] and N'Zi et al. [48]. BSDEs also appear in many fields of mathematical 
finance, see El Karoui et al. [27] or more recently El Karoui and Hamadene [26] for an overview. 
Problems such as pricing and hedging of European options (compare [27]), stochastic recursive 
utility (see [24]), utility maximization problems (e.g. [36]) and risk measures (e.g. [4, 50]) have 
been tackled using BSDE techniques. 

For linear BSDEs it is possible under certain intcgrability and boundedncss conditions to describe 
the first component of a solution as conditional expectation, compare [27] Section 2. This is already 
not possible anymore in the Lipschitz case and therefore solutions and their properties can mostly 
be studied numerically. In the meantime there is a huge literature on numerics for Lipschitz BSDEs, 
see [29, 10, 6] amongst many others. The present work focusses on BSDEs with drivers of quadratic 
growth which were first investigated by Kobylanski [41] in a Brownian setting and later extended 
to a continuous martingale setting by Morlais in [46]. Imposing certain growth and Lipschitz 
conditions on the generator and assuming bounded terminal conditions existence and uniqueness 
of quadratic growth BSDEs are guaranteed. However there exists much less research on numerics 
for quadratic BSDEs and only very few examples where an explicit solution is known. 

Motivated by this lack of examples and by the most important applications of this theory we 
analyze conditions under which one can find explicit solutions to a class of quadratic growth BSDEs. 
Our setting is as follows. We consider a forward affine process valued in the cone of positive 
semidefinite d x d matrices, and analyse BSDEs whose terminal condition and generator depend 
on this forward process. The main question from a BSDE-point of view addressed here is: Which 
structural conditions on the terminal condition and the generator arc needed (e.g. linear, affine, 
quadratic) to allow us to solve the BSDE explicitly? 
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We have chosen the forward process to be an affine process X on ST. As Kallsen and Muhle- 
Karbc [40] in the univariate case we can relate systems with Sj~-valued forward processes to mul- 
tivariate, realistically modelled utility optimization problems. Notice that affine processes have 
found a growing interest in the literature due to their analytic tractability which stems from the 
affine transform formula 



E 



e -Tt(X tU ) 



exp(-(t>(t,u) -Tr(iP(t,u)X )), 



for all t £ [0,T] and u £ Sf. The functions cf> and i(j solve a system of generalized Riccati ordinary 
differential equations (ODEs), which are specified via the model parameters. Affine processes have 
been applied in various fields in mathematical finance such as the theory of term structure of 
interest rates, option pricing in stochastic volatility models and credit risk, see e.g. [16, 22, 3, 23] 
and the references therein. Note that it is not necessary to look at matrix-valued affine processes 
as we do, but we could have equally chosen R" x R m , which was characterized by Duffie et al. 
[25], or R n x Sj~. For the sake of presentation we chose Sj~ since this state space is complicated 
enough to make important pitfalls visible (e.g. no infinite divisibility, no polyhedral property, etc, 
see [18]), but still allows for simple notation. We emphasize that most of our results carry over to 
the general state space case. 

The forward-backward system we consider consists of an affine process X on Sj~ and a BSDE 
whose terminal condition is an affine function of this process. Moreover the generator is allowed 
to have a more involved structure including a dependence on X and a quadratic dependence in 
the control process Z. In Theorem 3.3 we carry out which analytic form the generator and the 
terminal value of the BSDE need for the solution to be determined by a matrix ODE. This ODE is a 
generalized Riccati ODE which may explode in finite time because of its quadratic term. Therefore 
it is necessary to find conditions such that the ODE possesses a unique finite solution on the whole 
time interval [0,T]. 

We apply our results to the problem of maximizing expected utility of terminal wealth in multi- 
variate affine stochastic volatility models. This problem is typically approached either by stochastic 
control methods leading to Hamilton- Jacobi-Bellman equations or martingale methods which we 
will use here. The martingale method has also been used in [40], where the authors were able 
to solve the power utility optimization problem in several univariate affine models. They obtain 
the solution using semimartingale characteristics and represent the optimal strategies in terms of 
an opportunity process. Using a combination of martingale methods and our results on explicit 
solutions of BSDEs we extend these ideas to higher dimensions. In particular we derive explicit 
results for power and exponential utility in multivariate extensions of the model of Heston [34] 
and the model of Barndorff-Nielscn and Shephard [3] . The multivariate results presented here are 
mostly new and provide a thorough extension of Fonseca et al. [28] , where the power utility case 
in the multivariate Heston model is treated. We also want to mention the work of Leippold and 
Trojani [42] where in a multivariate setting the optimal strategy and the value function for power 
utility maximization is given. Their considerations are justified in our general affine setting. 

A particularly interesting application of our findings is the following: in the case of exponential 
utility we are able to provide analytic expressions for the indifference prices of variance swaps, 
which is well-known and can also be found in the literature. In the case of power utility - due to 
the additive structure of indifference prices one cannot find tractable expressions of those prices. 
However, the equally interesting concept of indifference value of change of numeraire is again 
analytically tractable. The indifference value of change of numeraire is the price one is willing to 
pay to swap one numeraire with another one. This can have two applications: one is the case 
where an institution actually bases their portfolio optimization, e.g., on fixed interest rates, even 
though interest rates are floating. The indifference value of change of numeraire is consequently the 
value of a swap contract particularly designed for compensating this model mispecification. The 
second one is a foreign exchange situation where replacing one numeraire by another one influences 
the optimal portfolio problem and therefore leads to an indifference value. We can provide fully 
tractable formulas in all these cases, see Section 4.3. 

In Section 2 we introduce necessary notation and collect several results by Cuchiero et al. [18] 
who give a complete characterization of affine processes on Sj~. Chapter 3 studies explicit solutions 
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to BSDEs and the connection with generalized matrix-valued Riccati ODEs. The results are then 
applied to expected utility maximization in Section 4. 

2. NOTATION AND CHARACTERIZATION OF AFFINE PROCESSES 

We start with the notation we use subsequently. We denote R+ = [0,oo), R_ = (— oo,0] and 

R++ = (0, oo), R = (—oo,0). We write C for the space of complex numbers. The space Md 

stands for d x d matrices with real entries and Id is the d x d-identity matrix. Sd is the space of 
symmetric d x d-matrices equipped with the scalar product Tr(xy). This inner product naturally 
induces a norm j| x \\= ^/Tr(xx). We write 5j" (or S~~) for the closed cone of symmetric d x d 
positive (or negative) semidcfinitc matrices and S^ + (or S d ) for the open cone of d x d positive 
(or negative) definite matrices. By dS d = S~£ \ S d + we denote the boundary of S d ■ The cones 
and + induce a strict partial order relation on Sd- We write x ^ y if y — x £ S d and x -< y 
if y — x £ Sj" . For i,j = 1, . . . , d, we also introduce the matrices e 4 -? £ M<j with efj = 1 and the 
remaining entries being 0. 

The Borel er-algebra on a space U C Sd is denoted by B(U) and frS'j' refers to the set of bounded 
real- valued measurable functions / on Sd- The vector space R d is equipped with the Euclidean 
norm | • |. We work on the finite time interval [0, T], where T > is fixed. 

Let us consider time-homogeneous Markov processes X with state space S d and semigroup 
(Pt)te[o,T] where 

Ptf(x)= [ f(0Pt(x,dt), 

and x £ / £ 65j" and p t a probability transition function. We refer to [51] for further details. 
The process X is not necessarily conservative. To construct a conservative process X we use the 
one-point compactification Sj" U A of <Sj~ . We correspondingly define 

p t (x, {A}) = 1 - Pt (x, S+), Pt (A, {A}) = 1, 

for all t £ [0,T] and x £ S d with the convention that functions / on S d are extended to Sj" U A 
by setting /(A) = and |A| = oo. Let us now define an affine process on St. 

Definition 2.1 ([18] Definition 2.1). A Markov process X is called affine if it satisfies the following 
conditions. 

(i) X is stochastically continuous, i.e. 

lim / f(Qp B (x,dQ= [ /(Oftfode), 
^Js+ J si 

for all / £ bSj and every t £ [0,T], x £ S#. 

(ii) The Laplace transform of X depends in an exponential affine way on the initial state. 
More precisely, there exist functions <j> : [0, T] X £j~ — ► R + and tp : [0, T] x S d — >• 5j~ such 
that 



p te -»(«0 = / e-^Mptfadt) = exp(-#t,u) - T*(V(*, «)*))> 
for all i e [0, T] and m, x £ Sj. 



In [18] Theorem 2.4 the authors establish the Feller property for every affine process X on S d ■ 
This permits a cadlag representation for a given process X. We may thus define the space fi of 
all cadlag paths w : [0,T] — > U A with oj(s) = A whenever u)(t— ) = A or = A for s > i, 
s,i £ [0,T]. For every x <G S^, F x is the (unique) probability measure on (O, \/te[o T] -^t^ ) sucn 
that P X (A = x) = 1, where (J^) is the natural filtration generated by X. We write J 7 ^ and 
respectively (J 7 ^) for the completion of Vte[o r] •^t' or (-^t^) with respect to P x . If we define 

T t = f| jf\ and J"= f) J-M, 

then the filtration (J 7 *) is right continuous and X is still a Markov process wrt. (J-*) for t £ [0,T]. 
In the following we will write "a.s." for "P x -a.s. for all x £ S d " . We call X a semimartingale if 
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X t l^ Xt( - s +^ is a semimartingale on (ft, J 7 , (J r t),^'x) f° r a U % G S^. For the definition of the charac- 
teristics of a semimartingale we refer to [37] Section II. 2. Note that semimartingale characteristics 
are always given wrt. a truncation function x '■ Sd — > Sd which is a continuous bounded function 
such that x(£) = £ for £ in a neighborhood of 0. 

Definition 2.2. We call (a, b, /3' J , m, /x, l, 7) an admissible parameter set associated to a truncation 
function x if it satisfies the following conditions. 

(i) The linear diffusion coefficient a belongs to the cone Sj". 

(ii) The constant drift term b is such that b y (d — l)a. 

(iii) The constant jump term m is a Borel measure on B(S~(f \ {0}) satisfying 

(|| £ || Al)m(d£) < 00. (2.1) 

'S+\{0} 

(iv) The linear jump term consists of a d x d-matrix /i = (/%•) of finite signed measures on 
B(S+ \ {0}) such that fx(E) G S+ for all E G 6(5+ \ {0}). The kernel 

M(M^«, (2 . 2) 

satisfies 

/ Tr(x(^)u)M(x, d£) < 00, for all i,ueS+ with Tr(xu) = 0. (2.3) 
■/sj\{0} 

(v) The linear drift coefficient is composed of a family (j8*-')i,j=i 1 ... 1 d of symmetric matrices 
with /3 i - J = /3 JJ G S'd for all i,j = l,...,d, and such that the linear map B : Sd —> Sd with 

d 

B(x) = J2 (2-4) 



fulfills 



Tr(B(x)u)- Tr( X (Ou)M(x,dO>0, (2.5) 

for all x, u G Sj" with Tr(xu) = 0. 

(vi) The constant killing rate coefficient t has values in R + . 

(vii) The linear killing rate coefficient 7 has values in S't. 

A full characterization of affine processes on S'j - is given in [18]. 

Theorem 2.3 ([18] Theorem 2.4). Let X be an affine process on Sf . Then there exists an admis- 
sible parameter set (a, b, f3 lJ , m, fi, 1, 7) wrt. a truncation function x such that the functions <j) and 
ip from Definition 2.1 (ii) solve the generalized Riccati ODE 

d(j>(t, u) 



dt 
dip(t, u) 
dt~ 



&(f/>(t,u)), 4>(0,u) = 0, (2.6) 

ay>(t,u)), 4>(o,u) = ues+, (2.7) 



with 



&(u) = Tr(6u) + L - (e-T'-tO _ l)m(df), 

J S+\{0} 

r e -Tr«) _ 1 , T r ( Y (C) u ) 

®(u) = -2uau + B*(u) +1 - H.II2 A1 MO, 

JS+\{0} II ? II A1 

where B*j(u) = Ti-(P lJ u) for i,j= l,...,d. 

Conversely, let (a, b, ft 13 , m, fi, l, 7) be an admissible parameter set associated to a truncation 
junction x- Then there exists a unique affine process on Sf and the condition of Definition 2.1 
(ii) holds for all (t,u) £ [0, T] x S* + , where 4> and ip are given by (2.6) and (2.7). 

Every conservative affine process on S* + with killing rate coefficients t = 7 = is a semimartin- 
gale. 
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Theorem 2.4 ([18] Theorem 2.6). Let X be a conservative affine process on 5j~ and (a, b, /3 IJ , to, fj,, 0, 
the related admissible parameter set associated to a truncation function X- Then X is a semimartin- 
gale whose characteristics (D, A, v) with respect to \ are given by 

D t = f (b+ f X(t)m(d0 + B(X S )) ds, 
Jo \ Js+\{o} J 

A t ,ijki = / Aijki(X s )ds, 
Jo 

v([0,t],G) = f (m(G) + M(X s ,G))ds, 
Jo 

for i,j,k,l £ {l,...,d}, t e [0,T] and G e B{Sj \ {0}). The matrix B is given by (2.4), M by 
(2.2) and A. l]kl by 

Aijki(x) = Xiuoiji + xuajk + XjkOLu + Xjia ik , (2.8) 

for all i,j,k 7 l € {1, . . . , d}, and x E Sf . Moreover there exists a dx d matrix of standard Brownian 
motions W such that X has the following canonical representation 

X t = x+ [ y r X~ s dW s 'Z+ [ E T dW s T V% ( 2 - 9 ) 




b + B(X s )+ / x{OM<%) ds 



X (0 (v x (ds,dO~Hds,dO) 

o Js+\{0} 

(H~x(0)v x (ds,dO, 

o Js+\{o} 

where £ € Md satisfies EE T = a and [i denotes the random measure associated with the jumps 
ofX. 

Remark 2.5. The canonical representation (2.9) follows via the canonical semimartingale represen- 
tation (see [37] Theorem II. 2. 34) and the construction of a matrix-valued Brownian motion. For 
the latter one has to find a matrix which coincides with the covariation of the affine process. 

Note that the constant drift term of an affine semimartingale is independent of the truncation 
function x while x influences the linear drift coefficient B. 

From now on we fix a truncation function %, and then write "admissible parameter set" for "ad- 
missible parameter set associated to truncation function x" ■ The affine process X with admissible 
parameter set (a, b, P tJ , m, /i, t, 7) is continuous if and only if m and pi vanish, i.e. (a, b, /3 y , 0, 0, t, 7). 
Since we only consider affine semimartingales we write (a, b, /3 IJ , m, [i) for (a, b, f3 13 , to, fj., 0, 0). 



To fix ideas let us give an example of a matrix-valued affine processes, the Wishart processes. 
These processes were first rigorously studied in [13] extending squares of matrix Ornstein-Uhlcnbcck 
processes. The dynamics of a Wishart process satisfies 

dX t = (b + HX t + X t H T )dt + ^/x~ t dW t Z + Y^dWj^fxl, (2.10) 

where &, H, S G Md and W is a d x d matrix Brownian motion. These processes have been widely 
used to model stochastic covariation, see e.g. [14], [20] and [21]. In order to obtain a well defnned 
matrix volatility process Bru required the constant drift part b = fcS T E for some k > d — 1. Then 
X has a Wishart distribution. In the above notation the admissible parameter set for the Wishart 
process is (E T E, b, /3 <J , 0, 0) with B(x) = Hx + xH T . 

In contrast to affine processes on the state space R™ x R™, which were fully characterized in 
[25], and where the diffusion term consists of a constant and linear part, the diffusion term of an 
affine process on with admissible parameter set (a, 6, , m, /1) only allows for a linear part of 



6 



ANJA RICHTER 



the specific form 

d 

u tj A tjk i(x)u u = 4Tr(xuau), x,ueSj. (2-H) 

Note that the necessity and sufficiency of conditions (m) and (v) in definition 2.2 was first shown in 
[18]. In particular formula (2.4) allows for a more general form than B(x) = Hx + xH T , x € , 
compare also [18], Chapter 2.1.2. 

3. Explicit solutions of quadratic FBSDEs 

In this section we examine how the solutions for a class of quadratic BSDEs can be reduced 
to solving ODEs. In contrast to many existence results in the literature, e.g. [49, 41, 45] and 
[5], where the generator / is usually required to satisfy certain Lipschitz and growth conditions, 
we suggest an analytic expression for / which gives the problem extra structure. Consider the 
following motivating example where such a form appears naturally. 

Take the one-dimensional Hcston model (compare [34]) for the dynamics of an asset H. The 
stochastic logarithm N of H satisfies 

dN t = rjRtdt + y/R~ t dQ u 

dR t = (b + \Ri)dt + (Js/R~tdWt, t e [0,T]. 

Here R is the stochastic volatility process, b, a > 0, r), A G R are constants and Q, W are two 
Brownian motions with constant correlation p £ [—1,1]. The volatility process R satisfies our 
definition of an affine process on R + with admissible parameter set ( \cr 2 , b, A, 0, 0). We study an 
investor who is interested in maximizing their expected utility from terminal wealth. The investor's 
initial capital is denoted byi€R and their trading strategies are deterministic functions tt of time, 
where ir(t) describes the amount of money invested in stock H at time t £ [0, T] . The wealth process 
X X ' 7T for initial endowment x and strategy 7r is given by 

X^ = x + f ^dH s = x + f ir(s)dN s , 

.In .In 



for t G [0,T]. We can solve the exponential utility maximization problem 

rX 

by finding the generator / of the BSDE 



V{x) = sup E [- exp {--fX^)] , 7 > 0, 



Y t = 0- I Z s dW s + I f(R s ,Z s )ds, (3.1) 



t Jt 

X.7T 



such that the process = — exp(— ^(X^ + Y t )), t 6 [0, T], is a supermartingale for all strategics 
7r and a martingale for a particular strategy n opt . The required generator can be shown (see Lemma 
4.13) to be 

f(r,z) = l(p 2 -l)z 2 + ^r 1 2 r-^pz^, (r,z)eR 2 , (3.2) 

similarly to [36] Theorem 7. Notice that the generator is quadratic in the z-component. In order 
to solve this BSDE we apply the Ito formula to an affine function of R. More precisely we make an 
affine ansatz for Y t = T(t)R t + w(i), t £ [0, T], where T, w : [0, T] — ► R are diffcrentiablc functions. 
This leads to 

T(t)R t + w{t) = T{T)R T + w(T) - J r(s)a y r R~ s dW s (3.3) 

T / , w x dT(s) dw(s)\ , 

T s 6 + XRs) + —^Rs + ds. 

ds ds J 

It can be immediately read off the equation that T and w must satisfy 

r(T) = w(T) = and Z s = T(s)a v ^R„ s e [0, T]. (3.4) 
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The finite variation parts of (3.1) and (3.3) coincide if for all s G [0, T] 

o = f(R s ,z s ) + r(s)(b + XR s ) + ^p-R a < dw{s) 



ds ds 

= l(p 2 - l)T 2 (s)a 2 R s + -Vi? s - -r 1P T(s)aR s + T(s)(b + XR S ) + ^R s + 
I Zj a 7 ds ds 

where we have used the equation for the generator / and formula (3.4) for Z. Equating coefficients 
this leads to an ODE of Riccati type 



dt 

with constants 



dr(<) - q r 2 {t) + ir{t) + c, r(T) = o, te[o,T], 



7 2/ 2 1 \ ; \ ^ „2 

-a (p -1), l = X apt], c=—^-q, 

2 7 



and an ODE of the simpler form 

dw(t) 
dT 

Hence the solution of (3.1) is 



T(t)b, w(T) = Q, te[0,T] 



Yi = T(t)Rt + J bT(s)ds 
Z t = T{t)a^fR t , te[0,T}. 



Generally, Riccati ODEs have the property that their solution can blow up in finite time, however 
our model parameter choices admit a non-explosive solution. In the one-dimensional case considered 
here we can even give an fully explicit solution, compare [12] Section 21.5.1.2. More specifically we 
distinguish two different cases depending on the value of 



1 \ 2 1 



d = I — Aqc = A — -ar]p + -2<7V(1 - P ) > 



7 J 7 2 



If d > 0, then 



Vd(T-t) _ i 

Tit) = -2c— =- — -, t e [o, ri. 



If d = 0, then p = 1, A = -crrj and hence 



7 

T(t) = ^ V 2 (T-t), te[0,T}. 
In both cases the martingale property of — exp(— r y(X x ' 7T P +Y)) then gives the value function and 



the optimal strategy 



V(x) = - exp ^-7 (^x + T(0)R o + bT(s)ds 



ir °*(j t ) = \ v -r(t)<T P , 



for x e R, t G [0,T]. 



In the previous example the ansatz Y t = T{t)Rt + w(t) and the method of equating coefficients 
enabled us to reduce the solution of a BSDE to solving ODEs. Now that we have seen how we exploit 
the affine structure in an one-dimensional example, we generalize this to BSDEs depending on affinc 
processes on S'j" and even an additional process which has affine semimartingale characteristics with 
respect to the affine process. The question is how general we are allowed to choose the generator 
and the terminal condition in order to still be able to apply the above method. 
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We associate the affine process X with admissible parameter set {a, b, , m, ^) to a BSDE. To 
allow for a more flexible financial modeling, especially in view of pricing of variance swaps, the 
BSDE will moreover depend on the matrix-valued process 

dO t =<j(t)^yx~ t dQt + {oi{t) + 02(t)X t )dt, te[0,T], (3.5) 

where 01,02 : [0, T] — > Md and a : [0, T] — > Md are continuous functions. Here the process Q denotes 
a d x d matrix- valued Brownian motion. This process will enable us to calculate indifference prices 
and delta hedges for variance swaps, see e.g. Section 4.3.2. 

Our (nonstandard) real-valued BSDE will have the following form 

Y t =F(X T ,0 T )- J Tr(ZjdW s )-J Tr(ZjdQ s ) (3.6) 

- I I K s (0((i x (ds,dO-v(ds,di;)) 
Jo Js+\{o} 

+ J f(s,X s ,Y s ,Z s ,Z s ,K s )ds, 

for t £ [0,T], where the terminal condition F is allowed to depend on the affine process X and on 
the process O. Recall that W is the Brownian motion of the underlying affine process X and the 
generator is a deterministic Borel measurable function / : [0, T] x 5 f | x R x Md x Md x M — > R. 

Definition 3.1. A solution to BSDE (3.6) is a family of processes (Y, Z, Z, K) £ K x M d x M d x R 
such that: 

(i) The equation (3.6) is a.s. satisfied. 

(ii) The stochastic integrals Z T ■ W , Z T ■ Q are well-defined. 

(iii) The integrability condition J Q T \f(s, X s , Y s , Z s , Z s , K s )\ds < 00 holds true. 

(iv) The mapping 1 i-> Yt is cadlag. 

(v) The process K satisfies !q J s +\{q} \ K t(0\( m (dO + M(X t ,d£))dt < 00. 

For a certain class of generators and terminal conditions we can give the solution to the above 
BSDE in terms of matrix ODEs. Suppose the terminal condition F : Sf x Md — >• M is affine, more 
precisely 

F(x,o) = Ti(ux) + Tr(ao) +v, x G Sj , o £ M d , (3.7) 

where u £ Sd, a £ Md and v £ R. Let us define the set L° as the space which contains all 
functions k : — > R. The class of generators / is more involved, more precisely the generator 
/ : [0, T] x S+ x K x M d x M d x L° -> R is allowed to have the following form 

f(t,x,y,z,z,k) (3.8) 

= Tr(zc zz (t)z T ) + Tr(zc zV ^(t)y/E) + Tr(c x (t)x) + c y (t)y + c t (t) 

+ Tr(zc zz (t)z T ) + Tr(zc zz (t)z T ) + Tr(zc 2 ^(t)y/x) 

g M (t,k{t)W{x,dt) 

si\{o} 

f 

sjAfo} 

for all (t, x, y, z, z, k) £ [0, T] x S+ x R x M<j x M d x L°. In the above, 

c zz j C z-Jxi c x , c zz , c zz , C z ^£ . [0,T] > Md, 
ct,c y : [0,T]->R, 

are continuous functions and gu '■ [0, T] x L° — > R is an M(x, <i£)-integrable function, x £ Sf, 
which is continuous in time. Finally 

9 z ^fxi 9x> 9a, 9zz, 9 z% /x '• 

[0,T] x L° — > M d , 



T r (zg z ^(t,k(0)V^) +Tiixg4t,k(0)) + 9t(t,k(0) + y9y(t,m)m(d0 
zg zz (t, k(0)z T ) + Tr(zg Bz (t, k({;))z T ) + Tr(%^(i, k(£))y/x)) m(d£), 
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gt,g v : [0,T] xL°^l, 
are m(d£)-intcgrablc functions which arc also assumed continuous in time. 

Remark 3.2. We restrict ourselves to the case where / does not depend explicitly on O since the 
structure of / is already quite involved. If this were not the case we would derive further coupled 
ODEs. However the terminal value F depends affine on O. We could allow for jumps in O provided 
those jumps have affine characteristics in X , but do not for reasons of brevity. It is possible to only 
consider functional forms of time cr(-) rather than 
a(-)VXm (3.5). 

We can now give the main theorem which describes the explicit form of the solution processes 
(Y, Z, Z, K) in terms of the solution to a system of generalized Riccati equations. 

Theorem 3.3. Let X be an affine semimartingale on Sf associated to the admissible parameter 
set {a,b, ft 13 , m, /i) such that 

[ || £ || (m(df ) + M(x, df )) < oo, x e S+. (3.9) 

"'IICII>i 

Furthermore suppose that there exists a unique solution T(-,u) : [0,T] — > Sd to the generalized 
Riccati ODE 

-^f^-=e(t,T(t,u)), T(T,u)=u, (3.10) 

with 

9{t,u) = 4uY, T c zz (t)Eu + &(t)u + u£' T (t)+B*(u) + c <f(t) (3.11) 
f Tr(^- X (0))+ gM (t,Tr( M 0) 

4+\ {0} iupai mj 

+ / f«S T 3 ^,Tr«)) +gT(t,Tv(uO)Xu + ug v (t,Tr(uZ)) + ^(i, Tr«)) 

JS+\{Q} \ 

+ C T T (t)ag i i{t 1 Tr{uZ))a T o-(t) + a T {t)ag £z {t 1 Tr(ut))^u 

+ ui: T gJ z (t,Tr(uO)a T a(t) + a T (t)ag z ^(t,TY(uO?jrn(dO, 

for (t,u) £ [OjT] x Sd- The functions J£{t) and ffit) are given by 

1 

tf(t) = c x (t) + a(t) T ac zz (t)a T a(t) + a T (t)ac z ^(t) + ao 2 {t), 
for all t 6 [0, T\. Let w(-,u,v) : [0,T] — > R be the solution of 
dw(t, u, v) 



= ^ y (t) + c^m + a(t) T ac iz m 



dt 

with 



w(t,T(t,u),w(t,u,v)), w(T,u,v)=v, (3-12) 



tu(t, u, v) = c y (t)v + Ct{i) + Tr(aoi (t)) + Tr(ub) 

+ [ (Tr«) + g y (t, Tr(tf))v + g t (t, Tr«))) ro(dO, 

for (t,u,v) G [0,T] x Sd x R. Then the above BSDE has the unique solution 

'Y t = Tr(r(t, u)X t ) + Tr(aO t ) + w(t, u, v), 

Z t = 2^X~ t T{t,u)Y> T , 

Z t =VX~ t a T (t)a, 

K t (£) =Tr(r(t, U )C), 

for all£eS+, te [0,T]. 



(3.13) 
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Proof. Let T and w be the unique solutions of (3.10) and (3.12). We apply Ito's formula for 
semimartingales to the function [0,T] x 5j~ x M d 3 {t,x,o) H> Tr(T(t,u)x) + Tr(ao) + w(t,u,v). 
Using representation (2.9) and (3.5) this gives 

Tr(r(t, u)X t ) + Tr(aOt) + w(t, u, v) 

= Tr(r(T, u)X T ) + w(T, u, v) + Tr(aO T ) -J 2 Tr (XT(s, u)^fx~ s dW t 
-jf (rT(T(s,u)b) + Tt(T(s,u)B(X s )) + Tt fr(«,«) y + 

Tr(r( S , U ) X (0)(M A "(^,^) - Kda.dC)) 



Mo} 



x(0"»(de) ^ 



7' 



* J S+\{0} 
T 



t JS+\{0} 



Tr(r( S , U )(£-x(0)M A (dfl.dC) 



Tr 



ar(y 

9s 



dw(s, u, v) 
Si 



(Z.s 



Tr ( acr(s) y/X s dQ 



Tr (a(oi(s) + o 2 (s)X,)) ds, 
where we have used basic properties of the trace to derive 

Tr (T(t, u)(V%dWt£ + S T dW t T V / ^t)) = 2 Tr f sr (^ u)y/X~ t dW t 



Because of the admissibility conditions (2.1) and (2.3) and the integrability condition (3.9) on the 
measures m and M, we may write the above equation in the following form 



Tr(r(t, u)X t ) + Tr(aOt) + w(t, u, v) 

= Tr(r(T, u)X T ) + w{T 7 u, v) + Tr(aO T ) -J 2 Tr (jT(s, u)^fx~ s dW s 

Tr(r(s, u)b) + Tr(r(s, u)B(X s )) + Tr (t(s, u) [ x(0™(<%) ) ) ds 

Tr(r( S , U )C)( M x (d S ,dC)-K^,dO) 
Tr(r( S ,u)(C-x(0)K^,dO 



(3.14) 



sj"\{o} 



t JS+\{0} 
T 



t JS+\{0} 



Tr 



dr(s,- 

ds 



-X, 



dw(s, u, v) 
ds 



ds 



Tr ( aa(s) \J X s dQ 



Tr (a(o 1 (s) + o 2 (s)X s ))ds. 



Hence the BSDE (3.6) is solved by (3.13) provided the finite variation parts of (3.14) and the BSDE 
coincide, i.e. if 



f(t,X t ,Y t ,Z t ,Z,Kt) 



s:\ioy 



-Tr 



Tr (r(t, - x(0)) + «)) 

s d + \{0} 

f dT(t, u) \ dw(t,u,v) 



\ dt 



-X t - 



dt 



Tr(ao 1 (t))-Tv{ao 2 (t)X t ), 
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for t G [0, T]. Using the special form (3.8) of the generator / and formulae (3.13) we calculate 
f(t,X t ,Y t ,Z t ,Z t ,K t ) 

= 4TY(^T t r(t,u)Z T c zz (t)i;r(t,u)^T t ) + 2Tv{^X~ t T{t,u)Y. T c z ^{t)^Y t ) + Tr(c x (t)X t ) 
+ Cy(t)Tr(T(t,u)X t ) +Cy(t)w(t,u,v) + c t (t) + Tv(y/Y t a T (t)ac zz (t)a T a(t)y/X~ t ) 
+ Tr( ^X~ t <J T (t)ac zz (t)2^T(t, u)y/x t ) + Tr( ^/x~ t cr T {t)ac z ^{t)^fx t ) 

g M (t,Tr(T(t,u)0)M(X t ,dO 

s:\io} 



SJ\{0} \ 



( Tr(2X t T(t, u)Y^g z ^{t, Tr(r(t, «)£))) + Tr(g x {t, Tr(T(t, u)fl)X t ) 

+ Tr(g y (t, Tr(r(i, u)0)r(t, u)X t ) + g y (t, Tr(T(t, u)0M*, u, v) 
+ g t {t)+Tr(^ t a T {t)ag zz (t,Tiir(t 7 u)0)a T a{t)^X~ t ) 
+ Ti(y/X~ t a T (t)ag zz {t, Tr(r(t, u)0)2Sr(t J u)y/X~ t ) 

+ Tr( ^X~ t a T (t)ag z ^(t, Tr(T(t, u)0) V^t)J H<%) 
- Tr(r(t, u)6) - Tr(r(t, u)B(X t )) 



S d +\{0} 



■&(T(t,u) X (0)m(de)- 



Tr 



I at 



dw(t, u, v) 
dt 



Tr(T(t,u)^- X (0))(m(dO+M(X t ,dO) 

si\{o} 

Tr(a 0l (t)) -Ti(ao 2 {t)X t )), 



where the last equality is obtained from (3.10) and (3.12), the definition of the adjoint operator 
B* and basic properties of the trace. □ 

Remark 3.4. Obviously our results also apply to 'standard' FBSDEs, where the BSDE is only 
allowed to depend the affine process X itself. 

In the remaining section we study the existence of generalized Riccati equations of the form 
(3.10). We always assume that (a, b, /3 U , m, n) is an admissible parameter set and that (3.9) holds. 
Furthermore we make use of the fact that the ODE (3.10) is equivalent to 

dT{t,u) 



dt 



6(t,T(t,u)), T(0,u)=u, 



(3.15) 



for (t, u) G [0,T] x Sd- This is derived via time reversal. The equations (3.10) and (3.12) are 
not coupled, in particular given T we need only to solve a linear ODE to get w. Thus, our main 
objective is to find sufficient conditions for existence and uniqueness of the solution T of (3.10). 

3.1. Solution of matrix Riccati ODEs in the continuous case. First we consider a simpler 
version of equation (3.10), more precisely 

dT(t,u) 



dt 



= 6{t,T(t,u)), T(T,u) = u, 



(3.16) 



8(t, u) = 4uE T c zz (i)Su + -Sf + u£? T (t) + B*{u) + c €(t). 

Note that for arbitrary u G Sd the solution of the above ODE may explode due to the quadratic 
term in the equation. In the corresponding continuous case with constant coefficients and terminal 
value u = 0we derive a closed form solution. We start with an assumption. 

(Al) Let the linear drift term B be of the form B(x) = xB + B T x, where B G Md- 

Proposition 3.5. Let (Al) hold and for t G [0,T] let the coefficients c zz (t) = c zz ,^(t) = 
Jtf, ^{t) = c (o be constant. Furthermore suppose a, c zz G Sf + and define the matrix-valued func- 
tions A : [0,T] — > M 2 d, Aij : [0,T] — > Md, i,j = 1,2, in the following way 



A(t) 



Au(t) A 12 (t) 
A 21 {t) A 22 {t) 



exp (T — t) 



+ B T -4£ T c zz £ 
-J£f - Bc x 
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Then, for the terminal value u = there exists a unique solution r(-,0) G S4 to (3.16) given by 

T(t,0) = A^(t)A 2 i(t), t€[0,T\. 

Proof. The result relies on theorems from Lie group theory and can be found in [33] Section 2.2. 
We only outline the procedure in which they linearize the matrix Riccati ODE. For every t G [0, T] 
let J(t) G M d be an invertible matrix and G(t) G M d . We set J(T) = I d , G(T) = and 

T(t,0) = J- 1 (t)G(t), t€[0,T], (3.17) 

which gives 

-d t G(t) = -d t (J(t)T(t,0)) = -d t (J(t))r(t,0) - J(t)d t (T(t,0)), 
and hence using (3.16) we have 

d t (J(t))T(t, 0) - d t G(t) =4G(i)£ T c 22 £L^, 0) + G(t)Jz? T + J(t)2T{t, 0) +J(t)BT(t, 0) + G(t)B T + J(t)tf. 
Equating coefficients, we get a system of 2d linear equations 

-d t (G( t ) j(t)) = (G(t) m(f T ^ BT ~*#_f)> ( 3 - is ) 

with terminal value (G(T), J(T)) — (0,ld). Via exponentiation we deduce its solution to be 
(G(t) J( t )) = (0 / d )exp((T-t)(^ T + ^ T -f^ C _f)). 

□ 

Remark 3.6. In the situation of Proposition 3.5 we can further simplify the computation of the 
component Y of the solution to BSDE (3.6) in Theorem 3.3. In order to determine Y it is typically 
not only necessary to solve (3.16), but also to numerically integrate Tr(T(t,u)b). This can be 
circumvented in many cases, compare [33] Chapter 4.2. Let Ct(t) = Ct, c y = and oi(t) = o\ for 
all t G [0,T]. With the notation from the above proof we have from (3.18) 



G(t) = (d t J(t) - J(t)J? - J(t)B) (4E 1 c zz E)- 
and hence, using (3.17) we get 

r(f, 0) = (j-\t)d t J(t) -££-B] (4£ T c 22 £)- 



Since 

w(t, 0, v) = v + J Tr(F(s,0)6)ds + J (c t + Tr(aoi)) ds, 

it follows with (3.13) for all t G [0,T] 

Y t = Tr(F(i, 0)Rt) + Tv{aO t ) +v- log || A 22 (t) || Tr ((4S T c 22 S)- 1 6; 

- (T - t) Tr f (J2f + B)(4E T c 22 S)- 1 6) +(T — t) (c t + Tr(aoi)) . 



To state further existence results we need additional assumptions, namely 
(A2+) c 22 G and^G S+ , 
(A2") c 22 G Sj" and ^eS^. 

Proposition 3.7. Lef (A2 + ) hold. Then, for every u € S'j", t/iere exists a unique solution T(-,u) G 
.Sj" to (3.16). Moreover, if 1 1— > (c 22 (t), Jzf^),^^)) is reaZ analytic in t G [0, T] and u G <Sj~ + , i/ien 

r(-,«)es++. 

Proof. For it G Sj" the existence of a solution r(-,w) G Sj" of (3.10) follows from [19] Theorem 
4.3. Indeed the continuous parameter set t i-» (-4E T c 22 (i)E, *«f(i)) with l{t)(u) = uJ? T (t) + 
Jif(t)u + B*(u), u G S'd, is admissible in the sense: -4Y, T c zz (t)Y,,tf(t) G Sj" and for all u G dSj", 
w G St such that Tr(uow) = and i G [0, T] we have 

Tr(l(t)(u )w) = Tt(u ^ t {t)w + ^(t)u w + B* {u )w) 

> Tr(u 3? T {t)w + £?(t)u w) 
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where we have used Ti(B*(x)y) = Tr (B(y)x) for all i,y£ Sd, as well as admissibility condition 
(2.5) and [18] Lemma 4.1. Since the right hand side of (3.16) is locally Lipschitz in u, uniqueness 
follows immediately. [19] Theorem 4.3 now implies the second assertion. □ 

Corollary 3.8. 7/(A2~) holds, then there exists a unique solution T(-,u) G S7 to (3.16) for every 
u G S7. Moreover, if t ^ (c Z2 (f), Jz? (t), ^(t)) is real analytic in t G [0, T] and u G S7~ , then 

Proof. Fix u G S7 . Finding a solution r(-,it) to the ODE (3.16) is equivalent to solving 

- ^ft^ =fo.r(t,-«)), ?(T,-u) = -u, (3.19) 
where T(t, — u) = —T(t,u) and 

(9(t, u) = -0(t, -u) = -4wE T c Z2 (i)Su + Jf(t)u + uJz? T (t) + B*{u) - *^(i) 

for (t,u) G [0,T] x S+. Assumption (A2+) and Proposition 3.7 imply that the ODE (3.19) has a 
unique solution T(-, —u) G 5j". Hence the unique solution of (3.16) is given by T(-, it) = — r(-, —it) G 
S7 . The second assertion follows immediately. □ 

3.2. Solution of generalized matrix Riccati ODEs. Before we derive sufficient conditions for 
the existence and uniqueness of equation (3.10) in the presence of jumps, we provide some properties 
of the function 9 given in (3.11). This enables us to use a comparison result for ODEs and prove 
existence of (3.10) with the help of the above existence results. We introduce quasi-monotone 
increasing functions in the sense of [54]. 

Definition 3.9. Let U C Sd be a an open set. We call a function 9 : U —> Sd quasi-monotone 
increasing if for all w,u G U and x G S'j" with Tr(iua;) = the inequality 

Tr((0(u> + u) - 9(u))x) > 

holds. Correspondingly we call 9 quasi- const ant if both 9 and —9 are quasi- monotone increasing. 

We work under the following conditions which ensure that 9, defined in (3.11) is quasi-monotone 
increasing. 

(A3 + ) For all t G [0,T], the functions gM(t,-), 9x{t,-), g z Jx(ti') an d 9y{ti') are non-decreasing 
on R + . For all (t,r) G [0,T] x R + , we have g M (t,r) G R+, g x (t,r) G Sj and 

S T 5.^(i, r) + gj^(t, r)S + g y (t, r) h 0. 

(A3 - ) For all t G [0, T], the functions gni(t, ■), g x (t, ■) are non-decreasing on R_ and the functions 
g z ^{t, ■), g v (t, •) are non-increasing on R_. For all (t, r) G [0, T] x R_, we have <7m(*, r ) & 
R_, g x (t,r) G SJ" and 

Z T g zV ^(t, r) + gj^(t, r)S + r) ^ 0. 

(A4 + ) For all i G [0, T], the functions g zz (t, ■), g Z z(t, ■), g z% /^(t, •) are non-decreasing on R + . For 
all G [0,T] x R + we have g zz (t,y) G 5j~ and 

^ T Wa.9^(i,2/)S + S T 3i,(i,y)a T ^(t) >= 0, 

°- T (*WiVi(*,j/) ^ o. 

(A4~) For all t G [0,T], the functions g zz (t,-), g z ^(t,-) are non-decreasing on M_ and g zz {t,-) 
are non-increasing on R_. For all (t,y) G [0, T] x M_ we have 

CT T (t)ag zz {t, y)S + £ T g zz {t, y)a T a(t) ± 0, 

o- 7 \t)ag Sy/ z(t,y) r< 0. 

Lemma 3.10. Suppose (A3 + ) and (A4 + ) /io/d. TTien i/ie map 9(t,-) is quasi-monotone increasing 
on Sj + for all t G [0, T] . 
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Proof. Consider an arbitrary t £ [0, T] and fix w, v £ + and r £ such that Tr(wr) = 0. First 
we consider the function 

and prove it is quasi-constant. Indeed, we have 

Tr (((w + w)4S T c„(t)E(tu + v) - «4£ T c 22 (t)£f) r) 

= Tr (w4£ T c zz (t)T,wr) + Tr (v4£ T c 22 0)£w) + Tr (wAY* 1 c zz (t)Y,vr) 

= 0, 

where the last equality follows from Tr(tur) = being equivalent to rw = wr = 0, compare [18] 
Lemma 4.1. With the same reasoning we have 

u^u^f T {t)+u^(t)+^(t), u£S+ + , 

is quasi-constant. Moreover due to the admissibility condition (2.5) the function 

Js+\{o> II 4 II Ai 

is quasi-monotone increasing. From the definition of fj, together with the fact that gu is non- 
decreasing in its spatial variable and Tr(xy) > for all x, y £ S'j", the map 

JS+\{0} II ? II Ai 

is quasi- monotone increasing. Due to (A3 + ) we have 

Tr ( r / ((w + v)Z T g zV ^(t, Tr((w + «)£)) + gj^(t, Tr((w + w)0)E(«> + «) 

- t;S T .g^(Tr(i,Tr«))) - ^(f, Tr«))£ V + («; + «)<?„(*, Tr({w + v)$) 

- vg v {t, Tr«)) + g x (t, Tr((w + «)£)) - Tr«)))m(d£)l 

= Tr(r / ( i ;E T ^(i,Tr(( W + W )0)+5 2 V(i,Tr(( W + W K))S« 
\ • / S+\{0} V 

- i;S T 3 2 ^(Tr(t, Tr«))) - gj^(t, Tr«))&, + vg y (t, Tr((w + «)£)) 

- vg y (t,Tr(vt;))+g x (t,Tr((w + v)t;)) - g x (t,Tr(v£)f)m(dZ)\ 

> o, 

which implies that 

Uh+ / (uS T ff2v ^(t,Tr«))+ ff T (t,Tr«))E M + ^((t,Tr«))+. g:c (t,Tr«)))m(^), 

is a quasi-monotone increasing function for all u £ . We finally obtain from (A4 + ) that the 
function 



u i — y 



(a T (t)ag zz (t,Tr(u£))a T a(t) + a T (t)ag Sz (t, Tr«))Su 

■S d +\{0} 



+«S T 5 T (i,Tr«))a T a(t) + a T (t)ag z ^(t,T r (u^)) m(d£), 

is quasi- monotone increasing for all u £ . □ 

In the following Lemma we show that the growth of 9(t, u) in u, where (t, u) £ [0, T] x S^, can 
be controlled under suitable conditions. 
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(A5 + ) The function c zz has values in Sj and gu and g x are of linear growth in the spatial 
variable, i.e. there exist continuous functions Cm '■ [0,T] — >• R++ and C x : [0,T] — !> 
such that 

g M (t,y)<C M (t)(\y\ + l), 
9x(t,y)lC x (t)(\y\ + l), 

for all (t,y) G [0,T] x R+. Also g z ^ and g y are bounded above in the spatial variable, 
i.e. there exists a matrix valued continuous function C z ^ : [0, T] — > such that 

9z ^{t,r) + g y (t,r) < C z ^(t), for all (t,r) g [0,T] x E+. 

(A5 _ ) The function c zz has values in and —31/ and —g x are of linear growth in the spatial 
variable, i.e. there exist continuous functions Cm '■ [0,T] — > R++ and C x : [0,T] — > + 
such that 

-9M(t,y)<C u (t)(\y\ + l), 
-g*(t t y) 1 C x (t)(\y\ + 1), 

for all (t,y) G [0,T] x R_. Also g z ^ and Pj, are bounded above in the spatial variable, 
i.e. there exists a matrix valued continuous function C z ^ : [0, T] — > such that 

g z ^(t,r)+g y (t,r) * C z ^(t), for all (t,r) G [0,T] x R_. 

(A6 + ) The functions g zz and <7 Zv /j, are of linear growth in the spatial variable, i.e. there exist 
continuous functions C zz : [0,T] — > such that 

9it{t,y) + g t ^(t,y) 1 C st (t)(\y\ + 1), 

for all (i, 1/) G [0, T] x R + . Also <? zz is bounded above in the spatial variable, i.e. there 
exists a continuous function C zz : [0, T] — > such that 

08*(*,lO r< Ci.(t), for all (t,y) G [0,T] x R+. 

(A6~) The functions — g zz and — <? Zv /^, are of linear growth in the spatial variable, i.e. there 
exists a continuous function C zz : [0, T] — > + such that 

-g zz (t,y) - g zVx -{t,y) ± C zz {t){\y\ + 1), 

for all (t,y) G [0,T] x R_. Also t/ zz is bounded above in the spatial variable, i.e. there 
exists a continuous function C zz : [0, T] — > such that 

ftta(t,lO r< C zz (t), for all (t,y) G [0,T] x R_. 

Lemma 3.11. Let (A5 + ) and (A6 + ) hold. Then there exists a continuous function K : [0, T] — > 
R ++ such that for (t,u) G [0,T] x Sj~ i/ie following inequality holds 

Tr(u9(t,u))<K(t)(\\u\\ 2 +l). 

Proof. For every (i, it) G [0,T] x 5j" we deduce that 

6»(£, it) ^ iiJz? T (t) + -Sfit + B* (it) + 1f(t) 

f C M (t) (Tr (<) + 1) + Tr - 
' / ii ^ 112 a 1 AH a sJ 

7s+\{o} II £ Ir ai 

+ / (iiS T C^(t) + Cj^mu + C x (t)(Tv «) + 1) 
Js+\{o} v 

+ a T (t)aC zz (Tr«) + l)a T cr(i) + cr T (»aC zz (t)£ii 
+ <r T (*)aC M (f)C& «) + l)Wd£), 



where we have used the Assumptions (A5 + ) and (A6 + ). In particular for it G St it then follows 
that 

Tr(u6(t, u)) < Ti(uu^f T (t)) + Tr(u_Sf + Tr(uB*(«)) + Tr(it^(i)) 
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m | / C M (t) (Tr( M fl + l)+Tr(u(g- X (fl)) 



Tr (u(uZ T C z ^(t) + CT(t)Zu + C x (t)(Tr «) + 1) 
5+\{o} V v 

+ a T (t)aC zz (Tr «) + l)a T a(t) + a T (t)aC Sz (t)T,u 
+ a T (t)aC zz (t)(Tr «) + 1)) )m(d£), 

which gives the result. □ 

The following theorem goes back to [54]. In the form stated here, it is adjusted to symmetric 
matrices. 

Theorem 3.12. Let U C Sd be an open set. Suppose 9 : [0, T) x U — > Sd is a continuous locally 
Lipschitz map with 9{t, •) quasi-monotone increasing on U for all t G [0,T). Given to G (0,T] let 
x,y : [0,io) U be differ entiable maps such that x(0) ^ y(0) and 

^ - 6(t, X (t)) < d A- 0(t, y(t)), t G [0, t }. 

Then we have x(t) ^ y(t) for all t G [0,to). 

We will use this theorem for proving an existence and uniqueness result for the generalized 
Riccati ODE (3.10). 

(A7) The mappings gM, 9z^fx-> 9xi 9yi 9zzi 9izi an d gz^/x are locally Lipschitz continuous in the 
spatial variable. 

Note that the map 9(t, •) may not be Lipschitz at the boundary dS^ , see [25] Example 9.3. So we 
look for a solution of (3.10) in St + . 

Theorem 3.13. Let the map t h-> (c zz (i), Jzf (<), t?(t)) be real analytic and suppose (A2 + ), (A3 + ), 
(A4+), (A5+), (A6+) and (A7) hold. Then the ODE (3.10) has a unique solution T{-,u) G S+ + 
for every terminal value u G S^ + . 

Proof. Since u h-> 9(t,u) is locally Lipschitz on + and 1 i— > 9(t,u) is continuous, standard ODE 
theory gives that there exists a unique local 5'j" + -valued solution T(t, u) of (3.15) for t G [0,i+(it)) 
with 

t + (u) = liminf {i > r II r(t,u) ||> n or T(t,it) G dSt] AT. 

Wanting a global solution we thus need to show that t+(u) = T. Since 9(t, •) may fail to be 
Lipschitz continuous at dS j~, we look at without its jump terms and define 

6{t, u) = AuY? c zz (t)T,u + u^ T (t) + 5?(t)u + B*(u) + <<f(t), 

for (t, u) G [0,T] x Sd- By assumption (A2 + ) and since the map t i-> (c zz (t), Jzf (t), ^{t)) is real 
analytic we may apply Proposition 3.7. Hence there exists a unique S*J" + -valued solution T of 

2EM = 0(t,f (*,«)), fl(o,«) = «, 

for all t G [0,T]. By means of (A3+) and (A4+) we have that for all (t,u) G [0,T] x Sj 
9(t,u) - 9(t,u) 

's:\{o} 



+\M II € II 2 Al 



(2 U £ T 5z ^(Tr(i, Tr«))) + 9x(t, Tr«)) + g y (t, Tr«))) m(d£) 

st\{0} 

(a T (t)ag zz (t, Tr«))a T a(i) + 2o T (t)ag zz (t, Tr«))£ U 
si\{0} 

+o T (t)ag z ^(t,Tr(uO))m(dO 
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to, 

and thus for all u G S d + and t £ [0, t+(u)) we obtain 

dT(t,u) - dt(t,u) - - 

— — 9(t,V(t,u)) h — 7^ 6{t,T(t,u)). 

Theorem 3.12 then yields for t £ [0,t+(it)) 

T(t,u)hf(t,u) e 

and hence t+(u) = lim inf n -+oc {t > : || T(t,u) ||> n} A T. By assumptions (A5 + ) and (A6 + ) we 
have from Lemma 3.11 that for all u £ S , j _+ and t £ [0, i+(u)) 

ft || T(t,u) || 2 = 2Tr(r(t, U )ftr(t,u)) < K{t) (|| r(i,«) || 2 +1) , 

for a continuous positive-valued function K. Then Gronwall's inequality applied to || T(t, u) || 2 +1 
gives 

II T(t,u) || 2 < e Jo K ^ ds (|| u || 2 +1) for t < t+(u), 

and thus t+{u) = T for u £ Sj + . □ 

Corollary 3.14. Let the map t (c zz (t), «£? (t), ^f(t)) be real analytic and suppose (A2~), (A3 - ), 
(A4 _ ), (A5 _ ), (A6 _ ) and (A7) hold. Then the ODE (3.10) has a unique solution T(-,u) £ S d 
for every terminal value u £ S d . 

Proof. Fix u £ . As we have already seen in the proof of Corollary 3.7, finding a solution 
r(-,w) to the ODE (3.10) is equivalent to solving 

_ arfo-tO = g (t f (*,-«)), f(r, -«) = -«, (3.20) 

eft 

where T(<, — u) = — r(i, u) and 
0(t,£t) = -6(t,-u) 

= -4uE T c^(i)Su + JSf(t)« + £LS? T (i) + B*(u) - tf(t) 
+ /■ Tt(fi^- X (0)) " 5A/(t, - Tr(^)) 
•^ d + \{o} 



e ii 2 ai 



•S d + \{0} 



+ o- T (t)ag 2z (t, - Tr(u^))Su - a T (t)ag £ ^(t, - Tr(u£)) )m(d£) 



for (t, u) £ [0, T] x S+. By Theorem 3.13 the ODE (3.20) has a unique solution f (-, -u) £ Sj and 
hence there exists a unique solution T(-, u) = — T(-, —it) £ S d . □ 

4. Application in multivariate affine stochastic volatility models 

In this chapter we apply the results of the previous chapter to the classical problem of utility 
maximization in a multivariate stochastic volatility setting. Stochastic volatility models are an 
extension of the Black-Scholes model, where the previously constant assumed volatility is now 
modeled as a stochastic process. The key feature of affine stochastic volatility models is that their 
Fourier-Laplace transform has an exponentially affine form. For a multivariate model consider the 
(i-dimensional logarithmic price process N whose stochastic volatility is given by an affine process 
R on S'j", then the following formula holds 



E 



e Tr(«_R t )+« 1 N t 



)Ra)+v T 7V +*(t 
e , 



for suitable arguments t £ [0,T], u £ Sd+iSd and v £ C d . The functions $ and ^ solve a system 
of generalized Riccati ODEs which are specified by the model parameters. This formula is the 
main reason for the analytic tractability of affine stochastic volatility models. In the multivariate 
stochastic volatility models mainly used in the literature, the dynamics of R follow 

dR t = {b + BR t + R t B T )dt + y/R~ t dW t Z + £ T dW T y/R~ t + dJ t , 
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Ro = reS+ 1 

where W is a matrix- valued Brownian motion possibly correlated with the Brownian motion driving 
N. Moreover b is a suitably chosen matrix in Sf , E, B are some invertible matrices and J is a pure 
jump process with a compensator that is affine in R. Without jumps this process is a Wishart- 
process (see also (2.10) and thereafter). They were introduced by [13] and have been applied to 
many different fields such as term structure modeling and derivative pricing in [31, 32, 20, 21]. 
The authors of [1, 2] consider multivariate stochastic volatility models for a class of matrix-valued 
Ornstein-Uhlenbeck processes driven by a Levy process of finite variation. Closer to our subject 
is the work of [28]. There the authors investigate the power utility maximization problem in a 
multivariate Hcston model where the covariation process follows a Wishart process. They obtain 
the optimal portfolio and utility via a duality approach. 

We start with the definition of multivariate affine stochastic volatility models. Thereafter we 
carry out the martingale property of the stochastic exponential of some process which will allow 
us to prove optimality in the utility maximization problem. More precisely we want to maximize 
expected utility of terminal wealth. The wealth process X x '™ is composed of the initial capital 
x 6 R and gains from trading with strategy 7r in the market. We want to solve the problem in 
presence of random revenues F which are paid at terminal time T, i.e. 

V(x) = sup E [U{X^ + F)] , x 6 E, 

where U is an exponential utility function. Once a notion of admissibility is fixed we call any tt G A 
an admissible (trading) strategy. Our aim is to explicitly describe the value function V and the 
corresponding optimal strategy 7r opt . Similarly we examine the problem 

V(x) = sup E [U(X*'*e F )] , 

where U now is a power utility function. It is known that the logarithmic utility maximization 
problem with F — can be solved explicitly for almost all semimartingalc models, see e.g. [30] and 
the references therein. This is why we do not consider logarithmic utility in this work. 



4.1. Definition. Multivariate affine stochastic volatility models are characterized via the joint 
Fourier-Laplace transform of the stochastic logarithm of the price process and the covariation 
process. We suppose that the discounted d-dimensional asset price process H is modeled as a 
stochastic exponential 

H t =H £(N) t , i€[0,T], 



where N is the discounted d-dimensional logarithmic price process with Nq = n E 1 
denote the stochastic covariation process with states in St and starting in Rq = r 6 Si 



Let R 



Definition 4.1 ([17] Definition 5.3.5.). We call a stochastic process (R,N) = (Rt, ^t)te[o,T] with 
values in Sj" x M. d a multivariate affine stochastic volatility model, if the following conditions are 
satisfied. 

(i) The pair of processes (R, N) is a stochastically continuous Markov process. 

(ii) Under the risk neutral measure Q, the Fourier-Laplace transform of (R, N) is exponen- 
tially affine in the initial states (r, n), i.e. there exist functions (t,u,v) <— > ^(t,u,v) and 
(t, it, v) i— > u, v) such that 



D Tr(uR t )+v T N t 



for all (t, u, v) G Q, where 



exp (Tr u, v)r) + v T n + $(t, u, v)) 



(4.1) 



Q= {(t,u,v) e [0,T] x S d + iS d 



,Tr(uR t )+v T N t 



< OO 



(iii) The asset price process H is a martingale under the risk-neutral probability measure Q. 
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4.2. The Martingale property. In the following sections it will play an important role under 
which conditions the stochastic exponential of a process involving an affine process becomes a true 
martingale. This problem has applications in fields including absolute continuity of distributions 
of stochastic processes (see [15]) and the verification of optimality in utility maximization as we 



use it here. On the state space '. 



and in a time-homogeneous setting the problem has been 



addressed already in [39, 44], which has then been extended to the state space S+ x R d in [17]. 

Suppose R is an affine process with admissible parameter set (a, b, /3 y , m, 0) associated with 
truncation function \ R - Let for all s G [0, T] 



{ITr^s)?)^!} 



and consider the process 



P t = / *J>{a)y/R a dQ a 



Tr 



'RAW. 



Tr ( (7q(s) v R s dQ, 



(4.2) 



(4.3) 



I / fe^C'MWO _ A d{n R (ds, dO - m{d£)ds),t G [0, T], 
o Js+\{o} v ' 

where ctq : [0,T] — > M. d and a\y , <7q , ct m : [0, T] M d are continuous functions of time. The 
d-dimensional Brownian motion Q is correlated to the matrix Brownian motion W by 



dQ t = dW t p+ s/l - p T pdD t . 

Here D is a d-dimensional Brownian motion independent of W and p a <i-dimensional vector with 
entries pi G [—1, 1], i = 1, . . . , d, satisfying p T p < 1. The process Q is another independent d x d- 
matrix Brownian motion. Finally p R denotes the random measure associated to the jumps of R. 
The goal of this section is to show that the stochastic exponential of P is a martingale which 
will help us proving optimality in the utility maximization problems considered in the following 
sections. 

Theorem 4.2. Assume (4.2), then the process £(P) is a martingale. 

It will be crucial that the conditional Fourier-Laplace transform of (R, P) with 

P = ]n(£(P)), 

is exponentially affine. This is stated in the following Lemma. 

Lemma 4.3. The conditional Fourier- Laplace transform of (R, P) has an exponentially affine 
form. More precisely, there exist functions (s,t,u,v) H- ^(s, t,u, v) and (s,t,u,v) H- $(,s,i,u, v) 
such that 



(Rs,Ps)] = cxp (Tr (V{s,t,u,v)R s ) + vP s + ${s,t,u,v)) , 



(4.4) 



for all (s, t, u, v) G X, where 

1 = {(s, t, u, v) G [0, T] x [0, T] x Sd + iS d x 

The functions <J> and 4" have the following form 
d$(s,t,u,v) 



s < t,E 



JT-c{uR t )+vP t 



(Rs,Ps) 



< OO 



ds 

d^(s,t,u,v) 



where 

&(s, u, v) = Tr{bu) + 



ds 



SJ\{0} 



&(s, *(s, t, u, v), v), t, u, v) = 0, 

M{s, *(s, t, u, v),v), *(t, t, u, v) = u, 



(4.5) 
(4.6) 



&(s,u,v) = 2uau + B*(u) + -v(v - 1) (2a Q (s)p T a w (s) + o^ v {s)a w {s) +o^{s)(Jq{s) 

+ cr g( s ) cr Q( s )) + vu{gq(s) P v + cr^(s))£ + vZ T (o- w (s) + pa^(s))u, 
for all (s, t, u, v) G T. 
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Proof. Since AP t > — 1 for all t <G [0, T] we have £{P) > and hence the process P is well defined. 
Ito's formula gives 



Pt = p t - \{p, P)t + M 1 + Ap *) - AP ° 



(4.7) 



Pt-^(P,P)t 



+ 



o Js+\{o} 
Tr(cr M (s)C) - e^*^) + l) m(d£)ds 



/0 ^S+\{0} 

Let us assume that there exist functions "J : X — > + iSd and <E> : Z — )• C such that the conditional 
Fourier-Laplace transform is of the form (4.4), i.e. ^(t, t, u,v) — u and $(t, t, u, v) = 0. We consider 
(s, i, u,v) £ I from now on and denote 

h(s,R s ,P s ) = cxp ^Tr(*(s,i,u, v)R s ) + vP s + <J?(s,i,u,u)J 

Then we obtain with Ito's formula for s < t G [0, T] 

dh(s,R s ,P s ) 

_ dh{ Sl R s ,P s ) 



s<t 



ids+j: 



dhjs^p) dh(s,R s ,p a ) j£) 

dr^ dp 



d 2 h{s,R s ,P s 



-d{Rij,Rki)s 
2 ■ irf i dr tj dr k i 2 



ld 2 h(s,R s ,P s 



dpdp 



■d(P,P) 



E 



d 2 h(s,R s ,P a ) 
dnjdp 



d{Rij,P'is 



j.STlh( R P \ h( f? p ^ V- dh(s,R s -,P s -) dh(s,R s _,P s _) 
+ / A h{s,R s ,P s ) - h(s,R s -,P s -) - } y Aifi,.., — AP< 



By the law of iterated expectations, for all s € [0,t] and s £ [s,t], 

h(s, R s , P s ) = E [e [e^^+^Kifc, P s )] P.)] = E [/i(s, ifc, A 

This means that /i(-,i?,P) is a martingale and hence the bounded variation part in the above 
equation has to be zero for all s < t E [0, T]. Straightforward computations then lead to 



= h{s,R s ,P s )[Tt 



d^(s, t, u, v) \ 9<E»(s, t, u, v) 



ds 



-R, 



ds 



2 Tr(P s *(s, t, u, u)a*(s, i, u, v)) (4.8) 



+ Tr((6 + B(R a ))V(s, t, u, v)) + -v{v - 1) Tr(P s (2a Q (s)p T o- w (s) + a^(s)a w (s) 
+ a^( S )a Q (s) +a Q (s)al(s))) + 2vTr(^(s,t,u,v)(a Q (s)p T +a w (s))E) 

+ [ fuTr(^(s)0 - we^Mfl + « + e Tr(*(«.*.«^)«+^(^WS) _ i) m (^) ds 
Js+\{o} v ' 

'- Tr(*(s, i, u, v)0 - v Tr(a M (s)0) m(^)ds ) . 



'S d +\{0} 

Here we have used the representation (2.9) of R, formula (4.7) and 

d(Rij,Rki) s = Aijki{R s ), 

d(P,P) s = Tr(R s (2a Q (s)p T a w (s) + a^(s)a w (s) + a^( S )<jQ( S ) + a Q (s)o-^(s)))ds, 

d(R,P) 3 = 2R s ((a Q (s)p T +a w (s))E)ds, s£ [0,T]. 

We obtain the ODEs (4.6) and (4.5) by equating coefficients. There exists a unique local solution 
for equation (4.6) since it is locally Lipschitz on S4 + iSd- From [52] Theorem 3.7 we know that 
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the time before the ODE possibly explodes coincides with the expectation on the left hand side of 
(4.4) to exist. This gives the result. □ 

In a time homogeneous setting [17] Theorem 5.3.4 gives necessary and sufficient conditions for 
the exponential of e p of a process P £ M. d to be a martingale. We adapt the proof to our time- 
inhomogeneous process P. 

Proof of Theorem 4-2. We show that e p is a martingale which gives the result due to e p = £(P). 
Since (4.2) gives 

e Ap >m(d£) = f e 1 *^^ m(d£) < oo, 

{|AA|>1} A|TrKW«)l>l> 

and P is a local martingale, [38] Lemma 3.1 implies that e p is a er-martingale. Due to e p > 
additionally we obtain by [38] Proposition 3.1 that e p is a supermartingale. Hence it is in particular 
integrable and thus (s, t, 0,1) £ X for all s < t £ [0, T]. From Lemma 4.3 we have for all s < t £ [0, T] 



E 



e p * 



{Rs,Ps 



exp (Tr(*(s, t, 0, 1)R„) + P s + *(a, t, 0, 1) 



Thus, ^(s, £, 0, 1) = and $(s, t, 0, 1) = for all s < t £ [0, T] would ensure that e p is a martingale. 
Indeed Lemma 4.3 implies 



&(s,0,l)= f fe- 1 ^««)-le^W«) + l-l)m(^) = 0, 
Jst\{o\ ^ ' 



's:\io} 

&{a, 0, 1) = B*(0) + ^1(1 - l)(2a Q (s)p T a w (s) + a^(s)a w (s) + a^a^s) + a Q (s)a^(s)) = 0, 

and hence the solution of 

d$(s,t, 0,1) 



ds 

a*(s,*,o,i) 



&(s,m(s,t,o, i),i) = o, ®(t,t,o,i) = o, 

^(s,*(s, t, 0, 1), 1) = 0, #(£,£,0,1) = 0, 



9s 

is *(s,£,0, 1) = and $(s,£,0, 1) = for all s < t £ [0,T]. □ 

4.3. Solution in a continuous multivariate afnne stochastic volatility model. In this sec- 
tion we introduce a continuous affine stochastic volatility model which is a natural multivariate 
extension of the Heston model. We then formulate the utility maximization problem in this model 
and solve it for power and exponential utility. This allows to describe the optimal strategy and 
maximal expected utility in terms of the model parameters and a Riccati ODE. 

Assume that there exists a financial market with one riskless bond with zero interest rate and d 
risky assets H = (Hi, . . . , Hd). The process H is modeled as stochastic exponential H = Hq£(N), 
where N = (Ni, . . . , No) is given by 

dN t = Rtrjdt + s/R t dQ u t £ [0,T], (4.9) 

with Q being a vector Brownian motion with values in M. d and 77 £ M. d being a constant vector. The 
stochastic volatility process is continuous and affine with admissible parameter set (a, 6, , 0, 0), 
i.e. it satisfies the SDE 

dR t = (b + B(R t ))dt+ ^/WtdWt^ + YJ dWj ^/Rt, R„ = r £ S+ , t £ [0,T], (4.10) 

where W is a d x d matrix-valued Brownian motion. The Brownian motion Q driving the assets 
returns and the Brownian motion W of the stochastic covariation matrix are allowed to be correlated 
in a certain way. Let p = (pi, . . . , pd) T be a vector with entries p; £ [— 1, 1], i = 1, . . . , d, and such 
that p T p < 1. With a d-dimensional Brownian motion D independent of W we can write 



dQ t = dW t p + y/l-p^pdDt. 

Hence the correlation between the scalar Brownian motions Q 1 and W mn is given by p n if i = m 
and else it is 0. The structure of the correlation between Q and W has been chosen in this way in 
order to ensure the model to be affine, which can be seen in the following Proposition. 
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Proposition 4.4. The pair of processes (i?, N) is a multivariate stochastic volatility model with 
functions $ and ^ solving 
d$(t,u,v) 



dt 

$(0, u,v) 
d^(t,u,v) 
dt 



6*(t, u,v), 
0. 

) + B*(V(t, it, v)) + (v - 7}) p T ^{t, u, v) 
+ *(t, u, v)Z T p (v - i]) T + -vv T , 



(4.11) 



(4.12) 

*(0,u,w) = u, 
for all (t, u, v) € Q. 

Proof. That (i?, AT) is stochastically continuous follows by definition. From Theorem 4.2 we obtain 
with o\v = o~q = (T M = and gq = —rj that the process £(X), where X t = — J Q 77 T \j R s dQ s , 
t G [0, T], is a martingale. This means that we can define a new probability measure Q by 
dQ = £ (X) T dP. The theorems of Girsanov and Levy then give that Q = Q + J Q ^/R^r/ds is a Q- 
Brownian motion. Hence N = J Q ^/R~ s dQ s is a local Q-martingale and again, according to Lemma 
3.1 and Proposition 3.1 in [38] the asset price process H = £ (TV) is a Q-martingale. 

It remains to prove formula (4.1). Since for all t e [0, T], the covariation d(Wij , X) t = 
—(y/Rtr])ipjdt, i,j = l,...,d,we have by Girsanov's theorem that 



W = W ■ 



r R s r\p ds, 



is a Brownian motion. Hence the dynamics of (4.9) and (4.10) under Q can be written as 
dN t = sfRtdQ u 

dR t = (b + B(R t ) - RtVP 1 ^ - Z T pri T R t )dt + y/WtdWt^ + YJ dWj ^f~R t . 



under the mea- 



Consider the conditional Fourier-Laplace transform 

sure Q for s < t. Since R and N are time-homogeneous processes, we have that for (t, u,v) E Q 
the process 



h(t-s,R s ,N s ) = 



^ uR ^ +vTNt \(R s ,N s ) 



is a martingale. We conjecture that the conditional Fourier-Laplace transform is of exponentially 
affinc form, more precisely that there exist functions 4 f : Q — >• Sd + iSd and $ : Q — > C d such that 
h has the form 

h(t - s, R S ,N S ) = cxp (Tr - s, u, v)R s ) + v T N s + $(< - s, u, v)) . 

Applying Ito's formula to hit — s,R s ,N s ) and taking into account that it is a martingale, the 
bounded variation term needs to satisfy 

/ d^(t- s,u,v) 



= Tr 



ds 



-R. 



\ d$(t—s,u,v) , , 1 -r 

j + g s + 2 T r(-R s *(* - s, u, v)a^{t - s, u, «)) + -v T R s v 



+ Tr ((6 + B(R S ) - R s T]p T Y, - S T prj T R S )V (t - s, u, v)) + 2 Tr(tf (t - s, u, v)R s Z T pv T ), 

for s < t. Here we have used the fact that d{Rij,N^) t = 2(p'S) i R t ji c dt for i,j,k = 1, . . 
Identifying coefficients leads to the following system of ODEs 

d^{t- 8,u,v) 



,d. 



ds 



= 2$(t - s,u,v)a^(t - s,u,v) + B*($(t - s,u,v)) 

+ (v — rj) p T S*(t - s,u,v) + \&(f - s,u,v)Y, T p(v - ?y) T + ^vv T , 



d$(t — s, u, v) 
ds 



M>(t- 8,u,v), 



with boundary conditions 



$>(t, u, v) = it, <f>(t,u, v)=0. 
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With a change of variable and since $ can be determined via integration as soon as \t is known, 
we derive (4.11), provided there exists a solution to the ODE for ^P, i.e. 



dt 

with 



6(u, v) = 2uau + B* (u) + (v — rj)p^Hu + uS T p(v — n) T + —vv T , 

for all (t,u,v) £ Q. Since 8 is locally Lipschitz, there exists a unique local solution of (4.12). 
From [52] Theorem 3.7 we know that the time before the ODE possibly explodes coincides with 
the expectation on the left hand side of (4.1) to exist. This gives the result. □ 

Remark 4.5. We have chosen to model the asset price process H as stochastic exponential of N. 
It is also possible to model H as ordinary exponential, i.e. H = H e N . In this case we have 
H = H £(N) with dN t = \fR~tdQ t + R t {v + \)dt, t £ [0,T]. Hence we are back in the setting 
considered above. 

4.3.1. Power utility. We set F = and assume that the investor's preferences are described by the 
power utility function 

U(x) = -X 7 , x > 0, 7 e (0,1). 

7 

Let A be the set of all <i-dimensional predictable processes tt that satisfy J Q T irjir s ds < oo a.s. Note 
that TTi denotes the fraction of the wealth invested in stock i, where i = 1, . . . , d. Any process it £ A 
is called an admissible (trading) strategy. Under these assumptions the wealth process evolves as 
follows 



Jo Jo Jo 

for t £ [0, T]. It can also be written as stochastic exponential 

X x t ^ = xE ttJ R s r]ds + J nJy^dQs^j , t £ [0, T]. 



As described earlier the investor wants to maximize their expected utility of terminal wealth. In 
order to model interest and exchange rates later, we want to take the function 

F(O t ) = Tr{aO T ), 

into account, where a is a d x d-matrix and Ot the final value of the process 

O t = f a^R~dQ s + f { 0l +o 2 R s )ds, t£[0,T}. (4.13) 
Jo Jo 

Here oi , 02 , a £ and Q is a d x d-dimensional Brownian motion independent of the Brownian 
motions W and Q. From now on we will write F a ' a 0l '° 2 for F(Ot) with Ot given by (4.13) to 
depict the structure of F in more detail. We want to solve the maximization problem 



V a ^' 0l ' 02 (x) = supE 



-(x*' 7r exp(F a << 7 ' 01 ' 02 ) 
7 v 



x > 0. (4.14) 



Our main result describes the value function and the optimal strategy explicitly in terms of the 
model parameters. 

Theorem 4.6. Let the linear diffusion term a belong to and suppose the linear drift term 

B in (4.10) is of the form B(r) = rB + B T r, r £ St , with B £ Md- Define the matrix-valued 
functions A : [0, T] — > Mid and Aij : [0, T] — > Md, i, j = 1, . . . , d, by 

A(t)= (An(t) A 12 (t)\ = ( ( t^W + ^t _ 2 a-^LETppT E \\ 

(4.15) 
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Then the value function and the corresponding optimal strategy are given by 

ya,a,o u o 2 ^ = exp ^Tr^ 1 (0) A 21 (0)r) + ^ Tr(A^(s)A 21 (s)b + ao 1 )dsj , 

< pt = j^- ^ + 2^ 22 1 (t)A 2 i(i)S T p) , x > 0, t G [0,T]. 

Before we prove the above theorem we motivate the present approach which can also be found 
in [36]. Note that in contrast to [36] the coefficients in the evolution of N arc not bounded. We 
solve this problem by using the martingale optimality principle, in particular we aim to construct 
processes as well as a strategy ir opt such that 

• L\ = U{X%* ex-p(F a > a ' 0l >° 2 )) for all n £ A, 

• L n is a supermartingale for all ir £ A and there is a particular strategy ir opt £ A such that 
L n ° pt is a martingale. 

Note that our assumptions on the filtration then imply that Lq = C for all ir £ A and a constant 
C > 0. Applying the utility function to I^' 1 exp(F a ' CT ' 0l <° 2 ) we get 

^(Xt* exp^' 7 ' 01 ' 02 )) 7 

= cxp uf jnjR^ds + jf jirJy/lT s dQ s - 1 jf 7 ^ T i? s 7T s d S + 7 F Q ^ 01 - ^ . 

This suggests the following choice of I/ 7r 

= x 1 cxp ^ ■yirjRsrjds + J jwJ^/R^dQ s ~\ j l n J R s ir s ds + Y t 

where Y is the first component of the solution of a BSDE with terminal condition ^F a ' <r ' 0l '° 2 , More 
precisely we want to find a generator / for the BSDE 

Y t = 7 F Q ^ 01 ' 02 - J Tr(ZjdW s ) - J Tr(ZjdQ s ) + j f(R s , Z„, Z s )ds, t £ [0, T], (4.16) 

such that its solution (Y, Z, Z) implies that meets the above requirements. 

Lemma 4.7. Let a £ St + , B(r) = rB + B T r with B £ M4 and recall (4.15). If the generator 
f : Sj x M d -^R is of the form 

f(r, M) = 5 Tr(zz T ) + 1 Tr(iz T ) + J _ + zp\\ (r, z, z) £ S+ x M d x M d , (4.17) 

then (4.16) is solved by 

Y t = Tr{A^{t)A 21 {t)R t ) + Tr(aO t ) + J Tr (A^ (s)A 21 {s)b + aoi) ds, (4.18) 

Z t = 2^R t A 22 1 {t)A 21 (t)Y, T , 

Z t = \f~Rto- 1 a, t£[0,T]. 
Moreover IF is a supermartingale for every strategy ir £ A and for 

< pt = 3^- (ri + 2A^(t)A 21 (t)^ T p) , t£[0,T], (4-19) 

the process L" P is a martingale. 
Proof. Let us define the constants 

Czz = \ Id + W L l) Ppr ' Cii = \ hu C ^=2(1^) WT ' Cx = W^y) W - 
Note that c zz is positive definite. By Proposition 3.5 we know that the ODE 

= r(t)S T c 2 ,Er(<) + B*{T{t)) + 2r(i)E T c^ + 2cJ^T(t) + ^a T aa T a + c x + ao 2 , 

r(T) = 0, 
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has the solution T(t) = A 22 l {t)A 2 i{t), t € [0,T]. We then obtain from Theorem 3.3 that the BSDE 
(4.16) with generator (4.17) is solved by (4.18). 

We show the local (super)martingale property for L" . With Ito's formula applied to L 7 * we have 
for all 7r e A 



LI (^TTj^RtdQt + Tr(ZjdW t ) + Tr(Z t T dQ t )) + ^ (j^jRtV ~ \l^t ' Rt*t ~ f(Rt, Z t , Z t ) 



dt 



1 ^Rt'Ktp Y + Z t + jy/l- p T (mJy/R t + Tv{Z t Zj))dt 



where we have used dQt — dWtp + y/T— p T pdD t . If the finite variation part satisfies dt <g> P-a.e. 



LI ( 77T, 1 R tV - -77r t ' Rt7T t - f(R t , Z u Z t ) + - 



y\/R t irtp + Z t 



y/l-pTpwJy/Rt 



< o, 



then we know that L^ is a local supermartingale. Indeed, since L™ > for all t <E [0, T] we only 
need to check whether 



1 



1 



-f(Rt,Z t , Z t ) < - 7 7r t ' RtV + gT** Rt*t 
This is equivalent to 

-f{Rt,Z u Z t ) < -^jRtV+^jirjRtnt - ± \Z t \ 2 



j\/Rt^tp T + Zt 



1 



y/l - p T pirj y/Rt 



1 



~fTr{y/R t ir t p T Zj) - -^irj R t ir t - -\Z t \ 2 



1 . 

27(1-7) 



Rtnt 



1 



1-7 



{y/RtV + Ztp) 



7 



2(1-7) 



'RtTj + Z t p 



\\^-\\z t \ 



If we use (4.17), we see that this inequality is true for every tt G A. For ir opt from (4.19) and 
applying the particular form of Z, the above inequality turns out to be an equality and hence the 
process L™ P is a local martingale. Note that ir opt S A. 

We proceed showing that is a true supermartingale for all tt G A. By definition there exists 
a sequence of stopping times (r n ) nS N converging to T such that L w At is a supermartingale. Since 
L" is bounded below by zero we may use Fatou's Lemma to pass to the limit: 



E[L?\F 8 ]=E 



lim LJ AT JF S 



< lim E [L? AT JF S ] < lim L? A = V s , s<te [0,T]. 



Note that from (4.18) and (4.19) we have 

Lf pt = x"<£ (fj (ir opt )Jy/R~ a dQ s + J Tr (2E(A 22 1 (s)A 21 (s)) T yfR~ s dW s ^ + J Tr (a T <JyfR~ s dQ t 

By choosing <tq(s) = 77r° pt , crw(s) = 2£(A 22 1 (.s)A 2 i(s)) T , ctq(s) = a T a and <r M = 0, s € [0,T], wc 
derive from Theorem 4.2 that L 7T ° pt is a true martingale. 

Proof of Theorem 4-6- Note that we derive from Lemma 4.7 for all ir G A 

E[U(X*'* exp(F a ^ 0l > 02 ))} = E 



□ 









"V 


E 


Lrp 


< E 






.7 . 




.7 °. 



1 



= -a; 7 exp(y )- 

7 

The strategy 7r opt is indeed optimal since we have that L^"^ is a martingale and hence 

T 



E 



U(X*>* exp(F a 



= E 



7 



This immediately gives the value function. 



□ 



Remark 4.8. In dimension d = 1, for the case F(Ot) = and with a slightly different choice of 
parameters, this result was derived by the authors of [40]. They represent the optimal strategy 
in terms of an opportunity process and use semimartingalc characteristics. In our setting the 
opportunity process is e Y , see also [47] and in particular [35] for a survey on the relationship 
between BSDEs and duality methods in utility maximization. On a heuristic level the result for 
d = 1 and F(Ot) = appears in [43]. Also using duality methods [28] derive a result similar to 
Theorem 4.6. 
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Finally we are able to give the indifference value of change of numeraire in two examples. Let 
us first look at the special situation where 

F -/ d ,o,03,o = _ Tr(o3)T or F-^ ' ol ^ = -/ Tr(o 1 +o 2 Rs)ds 7 o u . . . ,o 3 £ M d , 

Jo 

and understand this as the possibly stochastic discounting of the investors terminal wealth. The 
indifference value p of changing between those two numeraires is then defined by 

V- Id ' ' ol ' O2 (x-p(x)) = V- Id '°^°{x). 

Proposition 4.9. The indifference value of changing from a fixed interest rate F~ Id,0: ° 3: ° to the 
floating one F~ Id > > 01 ' 02 is 

p(x) = x - xexp ( - (tt{B^ (O)Bai (0)r) + jf Tr (B^ {s)B 2l {s)b - 03) ds 
-Tr(A£(0)A 2 i(p)r)- ^ Tr (A^ (s)A 21 (s)b - oi) ds jj , 

with 

^21 (i) ^22(«jy ^ V 2(T^)W °2 -lifW // 

B u (t) i?i 2 (t)^ ( (rr „ fe S W + ^ T -2a-^I] T W T E^ 



B 22 (t)J- ex »[( T *)^" 7 _ p X_ w T -^/S-B 
Proof. By Theorem 4.6 we have that 

y- Id '°' ol ' 02 (a;-p(a;)) = i(a; - p(^)) 7 cxp ^Tr^ 1 (0) A 21 (0)r) + jf Tr (^■ 2 1 («)A 3 i(«)6 - 01) ds^ 

V^ J - '°»'°(x) = exp ^Tr(B 22 1 (0)B 21 (0)r) + jf Tr (S 22 1 (s)B 2 i(s)6 - 03) ds^J , 

which gives the result. □ 

In a similar way we can describe the indifference value of change of numeraire from a fixed 
exchange rate 

F i d ,o,o 3 ,o = Tr ( 03 ) Ti 

with 03 € Md to a random valued exchange rate 

F a,a, 0l ,o 2 = Tr ( a o T ) j 

with a G and from (4.13). 

Proposition 4.10. The indifference value of changing from a fixed exchange rate F~ Id '°'° 3 ' to 
the floating one F a ' a > 0l >° 2 is 



xexp ^ ^Tr(S 22 1 (0)B 21 (0)r) + jf Tr (B^ 2 1 {s)B 21 {s)b + 03) 

-Tr(A 22 1 (0)A 21 (0)r) - / Tr (A^ 1 (s)A 21 (s)6 + aoi) ds j 



p(x) = x 



A 12 (t)\ _ / _ / t^W + -B T -2a - ^ pp^ 

B 2 i{t) B 22 {t)j y y -jhjvp s - b yy 

Since the proof is very similar to the proof of Proposition 4.9 we omit it here. 
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4.3.2. Exponential utility. In this section we want to solve the utility maximization problem for 
the exponential utility function 

U(x) = — exp(— jx), x G R, 

where 7 > denotes the risk aversion. Note that we have already discussed this problem in the 
one-dimensional case in the beginning of Chapter 3 and now study it in a multivariate setting in 
detail. We are also interested in pricing variance swaps which depend on the realized variance via 
utility indifference pricing. We consider the case where the variance swaps are not available in the 
market and the initial capital x is invested in the (incomplete) financial market H . For i = 1, . . . , d, 
a variance swap on the z-th asset of maturity T is a contract which pays 

1 f T 

- J (Ru) s ds 

at terminal time T in exchange for a previously fixed amount K{, That is to say the payoff of a 
variance swap on Hi is a function of Ot = J R s ds, more precisely 

F l {0 T ) =Tr(a ii O T )-K i , 

where a" = ^e %% . If we are only interested in the utility maximization problem without a random 
endowment, i.e. F % = 0, we define a" = 0, Ki = for i = 0. 

In this section we also need a notion of admissibility. For R™ x R"-valued affine stochastic 
volatility models Vierthauer [53] shows in Theorem 3.17 that the optimal strategy in the exponential 
utility maximization problem is a deterministic function of time. Motivated by this we introduce 
the set A of admissible trading strategies as the set of <i-dimcnsional deterministic functions of 
time 7T = (7r(£)) tg [ 0i T]- This time, the trading strategy ir describes the amount of money invested 
in the stocks H so that the number of shares is TTj/Hj for j = 1, . . . , d. The wealth process X x > n 
corresponding to strategy 7r and initial capital x is then given by 

X^ = x + Y] [ ^ldH %s = x + [ Tr T (s)R s r)ds + [ tt t ( S )^/R~ s dQ s . 

i=1 JO "i,s Jo JO 

Remark 4.11. Note that we measure the trading strategies n in different units than in the power 
utility case. This then leads to a similar exponential structure in the process L. This was also 
deployed in [36] for example. 

In the following Theorem we characterize the maximal expected utility from trading in the 
financial market in presence of a variance swap on the i-th asset 

V F \x)=supE[-cxp(-j(X^ + F 1 (O t )))], xeR, (4.20) 

ttEA 

and the optimal strategy ir F for i = 0, 1, . . . , d. 

Theorem 4.12. For i G {0, 1, . . . , d} let T l be the solution of the ODE 



dt 



T l (t) (-2 7 q + 2 7 E 1 pp 1 S) I«(f) + B*(T l (t)) (4.21) 

- -rm T PV T - w T sr(t) + -^ W T + a u , 

7 7 27^ 



r ; (T) = 0, 

for all t £ [0,T]. Then the value function has the form 

V F \x) = - exp ^-7 |x - lU + Tr(r J (0)r) + Tv{V \s)b)ds\ J , x e R, 

and the optimal strategy ir F is given by 

7r Fi (t) = ±r,-2r i (t)Y, T p, (4.22) 

for all t G [0, T] . 
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As before we want to use the martingale optimality principle in order to establish this Theorem. 
Therefore we dynamize the problem. For all i = 0, 1, . . . , d, we define 



■exp(-7(Xf •* + !?)), t&[0,T\, neA, 



where (Y 1 , Z % ) is the solution to 



Yj = F^Ot) - J* Tr((Zi) T dW s ) + ^ f(R s ,Zi)ds, te[0,T] 



(4.23) 



The generator / needs to be selected in a way such that L' K > % possesses the desired properties. This 
is done in the following lemma. 



Lemma 4.13. For i £ {0, 1, . . . , d}, let the generator f : x Md — > R have the form 

f(r, z*) = -I Tr(zV) T ) + ^ ~V^V - 7*V 
2 27 7 

Then the solution to (4.23) is given by 

Y t l = Tr{r(t)R t ) + Tr(a u O t ) - K { + J Tr(F (t)b)ds, 
Z\ = 2v^r l (t)E T , 



(4.24) 



(4.25) 



where T l £ is the solution to (4.21) and Ot = J^Rsds. Furthermore L n ' 1 is a supermartingale 
for every strategy n £ A and for 

(4.26) 



tt* (t) = - 2P(t)S 1 p, 



the process ' % is a martingale. 
Proof. Fix i £ {0, 1, . . . , d} and define 

Id), 



27 



2l 



3 VV 



Note that c zz is negative definite. Indeed, if p = 0, we have c zz = — ^Id £ S d . If p 7^ 0, we 
know that Id — pp T is the inverse of the positive definite matrix Id + 1 _ 1 p 1 p PP T and hence is itself 
positive definite. The conclusion is that c zz £ SJ~ . Then by Proposition 3.7 there exists a unique 
solution r l £ Sf. This allows us to find solution (4.25) via Theorem 3.3. 

Fix 7r £ A. By Ito's formula we see that IF^ can be described by the product of the local 
martingale 



Mr = -LY£ 1-7 
and the bounded variation process 



/ K T {s)^R s dQ s - f T r ((Zl) T dW s 
Jo Jo 



AT 



cxp 



(J (-lir T (s)R s v + lf(R s ,Zi) + ± 7 2 \VR s *(s)p T + (Z'l) T \ 2 



+ -j 2 \V^P T PK T (s)^ s \ 2 ds 



Theorem 4.2 implies that M 1 ' 1 is a true martingale. The process A*' 1 is non-increasing, if 

2 



1 ^{s)R sV + 1 f{R s ,Z i s ) + ^ 2 \^/R s K(s)p T + (^) T | 2 + l -^\^l^~p^ {s)^ > 



for all s £ [0, T]. This is equivalent to 



f(R t ,Zl)<^Tr(Zl(Zt) 



1 

7 



Rtn(t)-[ -py/Rtq-Zlp 



-^^V-jZlp 



2 \- T 
2 1 

2 



(4.27) 



±Tr(Zl(Zi) T ), 



EXPLICIT SOLUTIONS TO BSDES AND APPLICATIONS TO UTILITY MAXIMIZATION 



20 



which holds true by formula (4.24). Hence A*' 1 is non-increasing and L 71 "' 1 = M 7r 'M 7r '* is a su- 

pi pi- 

permartingalc. From (4.27) we see in particular that A* ' % = —1 is constant and thus L 7 * '* = 
— AP ' % is a true martingale. □ 

Proof of Theorem 4-12. Follows by the same reasoning as for Theorem 4.6. □ 

For all i G {1, . . . , d} , the indifference price of the variance swap 

F\O t ) =Tr(a"Or) - K 

on the i-th asset is defined as the value p l for which the investor is indifferent between buying F l 
for the amount p 1 and receiving a random income F l at terminal time T or not having it, i.e. 

V F \x-p l ) = V°(x), 

for all x G M. The optimal strategy tt f which attains the maximal expected utility in the presence 
of F l can be decomposed into a sum of a pure investment part ir° and a hedging component A 1 , 
i.e. 

ir F '(t) = ir°(t) + A l (t), t€[0,T}. (4.28) 
We therefore call A 1 the optimal hedge. 

Proposition 4.14. For i G {1, . . . , d} the indifference price p % and the optimal hedge A 1 of F 1 (Ot) 
are explicitly given by 

P l = -Ki + Tr((F(0) - r»)r) + f Tr((r*( S ) - r°(s))b)ds, 



A i (t) = 2(T i (t)-r°(m T p, t€[0,T], 
where T l and T° are the solutions of (4.21). 
Proof. Fix i G { 1 , . . . , d} and recall the value functions 

V F \x~p l ) = -exp f-7 {^x-p l - Ki + Tr(r l (0)r) + Tr(T l (s)b)ds 
V°(x) = - exp f -7 [x + Tr(r°(t)r) + / Tr(r°(s)6)ds J J , 



from Theorem 4.12. Equating them immediately gives the first part of the result. The second part 
then follows from (4.22) and (4.28). □ 

4.4. Solution in a multivariate affine stochastic volatility model with jumps. We now 

consider a model with jumps which is a natural multivariate extension of the model of [3] and has 
been applied e.g. in optimal portfolio selection, see [7] and the references therein. As before the 
asset price process H is modeled as stochastic exponential H = Ho£(N) with 

dN t = Rtrjdt + y/RtdQt, t G [0, T], (4.29) 

where Q is a d-dimensional vector Brownian motion and rj a constant parameter. By R we denote 
the Ornstein-Uhlenbeck-type stochastic process with dynamics 

dR t = (A + X(R t ))dt + dJ t , (4.30) 

and a starting value i?o = r. Here A £ S} and A : Sd — > Sd is the linear map A(r) = J2i j ^ r ij 
with /3 ij = /3 ji G S d and such that Tr(A(r)x) > for all r, x G S+ with Tr(rx) = 0. We denote its 
adjoint operator by A*. The process J is an independent affine process with admissible parameter 
set (0, b J , 0, m J , 0), starting at 0. Our goal is again to maximize the expected terminal wealth from 
trading in the market. 
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Proposition 4.15. The process (R, N) is a multivariate stochastic volatility model with functions 
$ and "J solving 



d$(t,u,v) 



dt 



= {\ + b J )V(t,u,v)- f ( e TV(£*(t,«,^)) _ l) m J (d£), $(0,u,«) = 0, (4.31) 

./.<7+\ Inl 



S^\{0} 



— = A*(*(t,u,v)) + -vv T , *(0,u,«) = u, 
/or a/Z (t, u, u) G Q. 

Proof. By construction (_R, AT) is again stochastically continuous. Moreover it follows that iJ is a 
martingale under the probability measure d<Q = £{— J Q r/ y/R 8 dQ s ) as in the proof of Proposition 
4.4. 

Again similarly to Proposition 4.4 we have that the conditional Fourier-Laplace transform for 
(t,u,v) G Q 

h(t - s,R s ,N s ) = E Q [e T ^ uR ^ +vTNt \(R s ,N s )\ , 

is a martingale. We assume that the conditional Fourier-Laplace transform is of exponentially 
affinc form, more precisely that there exist functions $:2^^ + iSd and $ : Q — > C d such that 
h has the form 

h(t — s,R s , N s ) = exp (Tr (V(t - s, u, v)R s ) + v T N s + $(t - s, u, v)) . 

We apply Ito's formula and stipulate that the bounded variation term needs to be zero. More 
precisely, this is equivalent to the equation 

„ / ty(t — s, u, v) „ \ <&(t — s.u, v) „ /, , , . „ . t , . nN 1 t „ 

= { ds V + ds + ^ + ( s) + ^ ~ s ' U ' w) ) + 2 sU 

( e ^(«*(*-«>«.-)) _ l)m J (^), 



'S+\{0} 

for s < t. Equating coefficients leads to the following system of ODEs 

\&(f — S,U,v) , 1 T . . 

— g — = A*(*(i - s, it, u)) + -i> u, ^(i, w, u) = u, 



* ,U,V) = (\ + b J )V(t~ s ,u,v)+ f ( e iv«*(t-,«*,«)) _ l) m J (^), =0. 



9s ^+\{o} 

With a change of variable we see that the above ODEs coincide with (4.31). The ODE for $ can 
be solved via integration, provided there exists a solution to 

d^it, u,v) 1 t 

y -^^=A*(^(t,u,v)) + -vv T , *(0,u,t;) = u, (t, u, v) G Q. 

This is a linear ODE, hence by e.g. [8], there exists a unique solution fy(t, u, v) G for all 
(t,u,v) G Q. This implies the result. □ 

AAA. Power utility. The investor's utility function is assumed to be 

U(x) = -x 7 , x > 0, 76 (0,1), 

7 

and we let F = 0. By „4 we denote the set of all d-dimensional predictable processes n that satisfy 
a.s. 7rJ ir s ds < oo. For i = 1, . . . , d, 7Tj again denotes the fraction of the wealth invested in stock 
i and any process 7r G *A is called an admissible (trading) strategy. Hence, for a trading strategy tt 
and initial capital x the wealth process has dynamics 

X?> n = x + f X*' w irJdN s =x + f X*' w irJ R s r)ds + f Xf *nj y r R~ s dQ Sl 
Jo Jo Jo 

for t G [0, T] which can also be written as a stochastic exponential 

Xf = xE (j wjR s r)ds + J irJ^R~ s dQ s ^J , i G [0, T]. 
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The investor wants to maximize their expected utility of terminal wealth, i.e. we search for the 
value function 

"1 



V(x) = sup E 



7 



, x > 0. 



We are able to describe the value function and the optimal strategy of the maximization problem 
in terms of an ODE. 

Theorem 4.16. Suppose the jump measure m J satisfies 

e -^ r ^m J (dO<oo, t€[0,T], 



'|Tr(r(i)£)|>l 

where T is the solution of the ODE 

dT(t) 



dt 

Then the value function is given by 



A*(r(*))- 



2(l- 7 ) 



W T , r(T) = o. 



(4.32) 



V(x) = ^x~> exp ( - Tr(r(0)r) - / 



Jo Js+\{0} 



for x > 0, and the optimal strategy n opt is 



opt 



1-7 



7], te[o,T] 



As before we solve the problem using the martingale optimality principle. Applying the utility 
function to X x '^ we get 

~(Xt' v y = ia; 7 exp [ [ jnj R s rjds + [ jnJy^dQs - \ ( 77rji? s 7r s ds J , 
7 7 \Jo Jo 2 Jo ' J 

for t G [0,T]. This suggests the following choice of L" 

LI = x 1 exp (^J jTrjR s rids + J ^tt] \f~R s dQ s - ^ J jirj R s ir s ds - F 4 J , 

for t G [0, T], and where Y is the first component of the solution of a BSDE with terminal condition 
0. More precisely we want to find a generator / for the BSDE 



Y t = - / / K s (£)((j, J (ds, d£) - m J {dH)ds) + / f(R s , K s )ds, 

Jt JS+\{0} Jt 

t G [0, T], such that with its solution (Y,K) L" satisfies the above requirements. 
Lemma 4.17. Let the jump measure m J satisfy 

f e- Tl '( r W«>m J (d£) <oo, te[0,T], 

where F is the solution of (4.32). Suppose the generator in (4.16) is of the following form 



(4.33) 



(4.34) 



/M) = ~2(l^) ??Tr7? 



si\{o} 



( e - fc «)-l + fe(0)m J (dO, (4-35) 



for all reS+andk: -> R. Tften SS'DS (4.33) is sofoed &y 



Y t = Tr(T(t)R t ) + / Tr(r(s)(6 J + A))da + 



t Js+\{o} 



(l- e -^( r WO) m J(^) rfs 



tf t (£) = ^(r(t)£), ie[o,T], £gS+. 

Moreover is a supermartingale for every strategy tt £ A and if 7r opt satisfies 



opt 



1-7 



r], te[o,T] 



(4.36) 



then L" is a martingale. 
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Proof. For all y £ R we define 



7 



2(l- 7 ) 



vv T , 9t(y) = -e y + l - y. 



Then we can see that by Corollary 3.8 there exists a unique solution T with values in Si to (4.32). 
As a result we find the above solution of (4.33) with Theorem 3.3. 
We apply Ito's formula which gives that for all n E A 

dLJ = VI ( 7 7r7 \/RtdQt) + V t (jnRtV - \Wt ' Rt*t + fifit, K t ) + ^j 2 irjR t n) dt 



LI 



SI\{0} 



st\{o] 



( e -K*«) - l) (fi J (dt,d^) - m J (d£)dt) 



l + K t (t))m J (<%)dt, 



where we have used integrability condition (4.34). This means that L* is a local supcrmartingale 
for all 7r € A, if the finite variation part dt (g> P-a.e. satisfies 

LI (^jRttj - ^jnj R t n t + f(R u K t ) + 'Rm 

+ I (e~ K ^~) _ i + K t (0) m J (dO) < 0. 
J5+\{0} V / ' J 

Since L" > 0, the generator / needs to fulfill 

f{RuK t ) < -jirjRti! + i 7 (l - "fWRm - f (e~ K ^ - 1 + Xi(0) m J (d£), 
which is equivalent to 

7 



f(Rt,K t ) < i 7 (i - 7 )|V^7r t - T^—.^Rtvl 2 - 



1-7 



2(1-7)' 



'RM 



s:\{o} 



--fft(e) 



With (4.35) this inequality is true for all 7r 6 A and hence L v a local supcrmartingale. Obviously 
the inequality is an equality for 7r opt , for which L 7! ° Fi is then a local martingale. 

Since L™ is bounded below by 0, the fact that L 11 is a supermartingale for all ir € A follows as 
in the proof of Lemma 4.7 by Fatou's Lemma. Note that the process 17 P is given by 



opt 



7 (n° pt ) T VRsdQ s + 



Js+\{o} 



( e -M0 _ l)( M J (d S) ^) _ m J {d£,)ds) , 



and hence, setting ctq(s) = j^zT), o--^(s) = and er^s) = — T(s), s 6 [0,? 1 ], we obtain the 

opt 

martingale property of L" by Theorem 4.2. □ 



Proof of Theorem J^.16. Note that we have from Lemma 4.17 

1 









"V 


E 


Lrp 


< E 


.7 L °. 




.7 . 







a; 7 exp(— Y ), n € A. 



Due to i 77 being a martingale we have 

S 



0' 



1 r 7r opt 

7 L ° 



and thus (4.19) is indeed the optimal strategy. This also provides the representation of the value 
function. □ 
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4.4.2. Exponential utility. We now show how to solve the utility maximization problem for expo- 
nential utility in the presence of jumps and random revenues F. The exponential utility function 
is given by 

U(x) = — exp(— 7X), 

where 7 > denotes the risk aversion. As before in Section 4.3.2 the deterministic d-dimensional 
functions tt = {^(t))te[o,T] form the set of admissible strategies A. For i = 1, tt 1 denotes 

again the amount of money invested in stock H l , where we recall that H = (H 1 , . . . , H d ) = £(N) 
with N defined in (4.29). In particular the wealth process X x ' n corresponding to a trading strategy 
7r and an initial capital x > satisfies 

" *(-) 



i=i 
ft 



X^=x + J2 / W 1 dH l , s =x + / M^dK 



[ ir T (s)R s r]ds+ [ TT T (s)y r R s dQ s 
Jo Jo 



As in Section 4.3.2, using the same notation, we will also compute the utility indifference prices for 
variance swaps. This means that we need to solve the problem 

V pz (x) = sup E [U + F*)] , x > 0, 

7TGA 

for i £ {0, 1, . . . , d}, which is done in the following Theorem. 
Theorem 4.18. Let i £ {0, 1, . . . , d}, 

[ e^ T ^ r '^m J (d^) < 00, se [0,T], 

./|-yTr(r*(a)£)|>l 

where T 1 is the solution of the ODE 

^=A*(r i (i)) + - W T +a», F(T) = 0. (4.37) 

T/ien £/ie waZue function satisfies 

V pi (x) = — cxp f — 7 fa; — A'j + Tr(r*(0)r) + / Tr(P(.s)(&' 7 + A))ds 



- ( T [ (er^ Ti ^ -l)m J (dOds 
1 Jo Js+\m v ' 



1 Jo JS+\{0} 
and the optimal strategy ir F is given by 

Tr Fi (t) = -r], t€[0,T]. 

1 

We will use the martingale optimality principle again and construct for every i £ {0, 1, . . . ,d}, 
a process 

L^ = -exp(-7(Xr+Y/)), t€[0,T), tt € A, 
where (Y l ,K l ) is the solution of 

V = F i - ( [ Ki(0(ti J (ds,dt)-m J (dZ)ds)+ [ /•:/,'-. A'l :■«/>. (4.38) 
Jt Js+\{o} Jt 

with a generator / such that L*'' 4 satisfies the conditions 

• the terminal condition L^' 1 = U(X^,' V + F l ) is satisfied for all tt £ A, 

• the process L"' ! is a supcrmartingale for all tt £ A and there is a tt £ A such that L v 
is a martingale. 

The following Lemma shows how the generator of BSDE (4.38) needs to be chosen in order to 
meet the above requirements. 



34 



ANJA RICHTER 



Lemma 4.19. Let i € {0, 1, . . . , d}, 

e ^ r< ««W(^) < oo, 



| 7 Tr(r>(t))|>l 

for all t £ [0, T] and with T l being the solution of (4.37). Let the generator f in (4.38) have the 
form 

f(r, fc l ) = ^- v T r V - - [ (e? fc '«> - 1 + 7 F(£)) m J (<%), (4.39) 
2 7 7 Js+\{o} v y 

for all r € iSj - and k l : Sf — > K. TTien t/ie solution of BSDE (4.38) is given by 



Yj = Tr{r l {t)R t ) + Tr (a u J R s ds^j -K 4 + Tt(T i (s)(b J + X))ds 
Lf*WQ _ m J (dt)ds, 



1 ' T 



(4.40) 



1 Jt Js+\{o} 

Moreover for all ir E A the process LF' % is a supermartingale and L" is a martingale. 
Proof. Fix i E {0, 1, . . . , d} and define for all y E R 

9x = ^-W1 T + a", g t (y) = -- (e™ - 1 + -yy) . 

Using Proposition 3.7 we see that there exists a unique solution T l E Sj" of (4.37). This implies 
(4.40) by Theorem 3.3. 

Fix 7r E A. Note that L"' 1 can be written as a product M n ' l V 7T ' t of the two processes 

M?> d = -L^st-j f K T {s)^R~ s dQ s + f [ - l) (v R (ds, d£) - m J (d0ds) , 

V J o J o Js+\{o} v ' ) 

Vr = - cxp (f (^- 1 7T T (s)R s r 1 + !f{R s ,K) + \-y\ T {s)R s -k s 

+ ^ o} (e^® - l + 7^(f)) ™ J (^)) ^ ■ 

Setting <jq(s) = — 77r s , aw{s) = and cr mu (s) = 7 r(s), we have from Theorem 4.2 that M 7r ' 1 is a 
true martingale. In order for V*' 1 to be decreasing, it needs to be ensured that 

- 1 Tr T (s)R s 7 1 + 1 f(R s ,Kl) + i 7 V( s )i?^ s + f (e 7 ^(?) - l + iKi(t)\ m J '(df) > 0, 

2 • y s+\{0} v 7 

(4.41) 

(is ® P-a.e. Taking formulas (4.39) and (4.40) into account this is indeed true, since (4.41) is 
equivalent to 

-f(RuKl) < i 7 kW T y/Rt- -V T VRt\ 2 - ^ T RtV 
z 7 z 7 



f I(e^«)-H- 7 JiC<(0)m J (de). 
•/stutrt 7 v y 



/s d +\{o} 7 

Since M"' 1 is a martingale and U 17 -* is non-increasing, L 71 ''* = M^^V^' 1 is a supermartingale. It is 
straightforward that V£ ' l = —1 for s £ [0, Tj and thus L 77 ,l = —M 71 - % is a true martingale. □ 

Proof of Theorem 4- 18. The proof follows the same reasoning as the proof of Theorem 4.16. □ 

Recall that for i € {1, . . . , d} the indifference price of the variance swap F l on the i-th asset is 
the value p l such that for all x G R the value V F (x — p 1 ) equals V°(a;). 
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Proposition 4.20. For i £ {1, . . . , d} the indifference price p l is explicitly given by 



T 



p' = -Ki + Tr((P(0) - r°(0))r) + / Tr((r (s) - T°(s))(b J + \))ds 

Jo 

e 7Tr(r-( S )0_ e 7TV(r°M?) N ) m ./ (de); 



1 ? T 



JO 

( e 7Tr(r>)£) _ 

7 Jo Js+\{o} 
where T 1 and T° are the respective solutions of 

^ii = A*(T*(t)) + — m T + a« r*(T) = 0, 

^=A*(r°(t)) + i- w T r°(T) = o. 

Proof. For i S {0, 1, . . . , d} it follows from Theorem 4.18 that the value functions have the form 
V F *{x-p i ) = -exp f-7 ( 2 ;-p 4 -K l +Tr(r(0)r)+ / Tr(F(s)(6 J + A))da 



n 



1 

1 Jo Js+\{o} 



Equating y F '(a; — p J ) and y°(a;) for i = 1, . . . , d, immediately gives the result. □ 
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